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PROGRESS REPORT 


I. Publications 

The following papers have been completed during the past six months 
and represent the progress made in describing and interpreting coronal 
plasma processes and the relationship between the Solar corona and its 
stellar counterparts: 

"On Stellar X-ray Emission", by R. Rosner, L. Golub and G.S. Vaiana, 
is a review article which will appear in Annual Reviews of Astronomy 
and Astrophysics . This article describes the manner in which Solar 
Physics, Stellar Physics and General Astrophysics have joined together 
into the entirely new discipline of stellar x-ray astronomy. A central 
issue throughout the review is the extent to which extrapolations from 
solar physics are justifiable and the manner in which stellar observa- 
tions are being used to augment and test the solar studies. Topics 
discussed include: "X-ray Emission From Stars", "The Solar Analogy and 
Modelling of X-ray Emission", "The Physics of X-Ray Emission from Solar- 
Type Stars", and "The More General Astrophysical Context". The paper 
concludes by outlining the most pressing problems in coronal physics 
and discussing prospects for their solution. 

The paper "EINSTEIN X-ray Survey of the Pleiades", by G. Micela 
et al . has been accepted for publication in the Ap. J . This paper rep- 
resents a significant step in understanding the relationship between 
the interior structure of solar-type stars, the presence of dy lamo acti- 
vity and the type and strength of coronal activity which results. Prior 



to these observations, a clear correlation had been observed between 


the strength of x-ray emission from solar-type stars and the rotation 
rates of the stars. The link has been presumed to be via strong inter- 
nal differential rotation which causes vigorous magnetic field produc- 
tion and results in the obscirved x-ray activity. Young stars in general 
rotate rapidly and have high levels of x-ray emission; the only excep- 
tion has been stars earlier than A7, which do not have significant outer 
convection zones and are not detected in x-rays. 

The Pleiades observations show an apparent contradiction with this 
picture, because although the G-stars in the cluster are strong x-ray 
emitters, the K-stars are relatively weak in spite of their very rapid 
rotation rates; moreover, few if any M-dwarfs are detected even though 
most nearby M-dwarfs are strong and highly variable x-ray sources. 
The key to the puzzle appears to be the fact that the Pleiades is a 
very young cluster and that the K and M stars may not yet be evolved 
onto the main sequence. The G-dwarfs are just old enough to have set- 
tled into a structure with an outer convection zone and a rapidjy rota- 
ting core. This interpretation strengthens the case for those solar 
dynamo models which involve magnetic field amplification at the base 
of the HCZ, at the shear interface with a radiative core. 

The paper "Closed Coronal Structures VI...", described in our last 
progress report was revised and has appeared in the Ap. J . In addition, 
the paper, "X-ray Survey of Main-Sequence Stars With Shallow Convection 
Zones" has appeared in the Ap. J . and discusses the onset of convection 
in solar-type stars and its relation to the dynamo and magnetic field- 


related activity. 



The paper "Implications of the 1400 MHz Flare Emission..." by Holman, 

Bookbinder and Golub, was presented at the conference on Stellar Con- 

Tinuum Radio Astronomy in Boulder, CO. The paper compares solar and 

stellar radio flares and discusses the possible emission mechanisms by 

13 

which the high brightness temperature of > 10 K could be produced 

by solar-type flares. Predictions of emission properties such as polari- 
zation, variability timescales and source size can be coupled with measure- 
ments of x-ray emission and magnetic field strength in order to test 
Solar flare models in a more extreme parameter regime. 

A paper in progress, "On Magnetic Field Stochasticity. . ." by Antonucci, 
Rosner and Tsinganos uses SMM data on line broadening during Solar flares 
in order to test a radical new theory of Solar flare triggering and develop- 
ment. Results indicate that flares may be caused by a process which 
leaves the external appearance of loops nearby unchanged (initially) 
but produces completely disordered field line structures within the loop. 
Microscopic turbulence and globally incoherent heating would then be 
properties of the initial flare onset. This view seems to explain many 
of the previously puzzling features of the BCS and SXP data. 


II. Talks, Abstracts and Invited Presentations 

"An X-ray Survey of Solar-Type Stars" by J. Bookbinder, P. Majer, 
L. Golub and R. Rosner, 165th AAS Meeting, Tucson. 

"Correlation Between X-ray Luminosity and Rossby Number", by G. 
Mlcela, S. Sciortlno and S. Serio, Symposium "X-ray Astronomy In 1984", 
Bologna, Italy. 

"Optimized Variability Analysis in Non-Periodic Sources", by A. 
Collura, S. Sciortlno, A. Maggio, S. Serio and R. Rosner, Bologna Sympo- 
sium. 

"X-ray Variability in K and M Dwarfs", by C. Ambruster, S. Sciortlno 
and L. Golub, abstract to 166th AAS Meeting. 

"The Time Resolution Domain of Stellar Radio Astronomy", by J. 
Bookbinder, paper to appear. 

"Theory Tested by Means of the Stars", by L. Golub, NASA CP-2358, 
eds. Underhill and Michalitsianos . 

"Implications of the 1400 MHz Flare Emission..." by G.D. Holman, 
J. Bookbinder and L. Golub, in Radio Stars, eds. Gibson and Hjellming, 


Proc . Conf. on Stellar Continuum Radio Astron., Boulder (1984). 
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ON STELLAR X-RAY EMISSION 


R. Rosner 

Department of Astronomy, Harvard University, Cambridge, MA 02138 
L. Golub 

Harvard-Smithsonian Center for Astrophysics, Cambridge, MA 02138 
C. S. Valana 

Osservatorio Astronomico, Kalermo, Italy 


I. INTRODUCTION 


The past five years have seen the emergence of an entirely new 
astronomical discipline, stellar x-ray astronomy. It lies at the crossroads 
of a number of apparently disparate astronomical research areas: 

(a) Solar physics : Stellar x-ray astronomy constitutes the major testing 
ground for theories of stellar surface activity which have been developed over 
the past two decades in the solar context. Stellar observations allow one to 
test predictions regarding the dependence of activity level on stellar 
parameters which cannot be varied in the course of observations restricted to 
the Sun. These parameters Include stellar age, the mean stellar rotation 
rate, the depth of the convection zone, and stellar surface characteristics 
such as surface effective temperature, effective gravity, and elemental 
abundances . 

(b) Ste^^ar p h ysi cs: The stellar surface activity which underlies the 
observed x-ray emission provides the means by which stellar splndown and mass 
loss occur. Until the advent of x-ray observations, other proxy indicators of 
stellar activity (such as Ca II H and K emission: Wilson 1966, 1978) proved to 
be insufficiently sensitive to provide the large volume-limited samples of 
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stars throughout the Hertzsprung-Russell diagram which are required for 
systematic and unbiased study of the full range of possible stellar activity 
levels (as opposed to study of only the most active examples). The presence 
of such surface activity also has Important implications for more traditional 
studies of stellar properties, such as determinations of surface elemental 
abundances and stellar age (Duncan 1981; Soderblom 1983, 1984); this occurs 
because surface activity on the Sun is now Icnown to modulate the measurement 
of abundances for important species such as Lithium (cf., Giampapa 1983). 

(c) General Astrophysics: Because of the particle and mass ejections 
which accompany stellar surface activity, stars provide one of the main 
sources of mass input to the Interstellar medium (ISM), and may provide at 
least some of the input source particle spectrum for galactic cosmic rays (see 
review by Montmerle 1984). The observed x-ray emission not only provides a 
primary diagnostic for the presence and level of such activity (and in 
principle allows determination of the elemental abunoances of the source 
material for the ejected matter), but also figures as an important source of 
ionizing flux in the Immediate stellar neighborhood (especially in regions of 
star form'itlon). It is also thought that stellar activity levels modulate the 
accretion rate in low-mass compact binary systems such as cataclysmic 
variables (Rappaport, Verbunt, and Joss 1983; Spruit and Ritter 1983: 
Patterson 1984); that the x-ray emission from young stars may play an 
Important role in the chemical state of the ISM in star-forming regions (viz., 
Feigelson 1983); that the x-ray emission from low-mass main sequence stars can 
contribute significantly to the galactic component of the diffuse soft x-ray 
bac)<ground (Rosner et al . 1981; Boo)<binder et al . 1984): and, finally, that 
the physical processes which figure in stellar surface activity (e.g.. 
magnetic dynamo activity, magnetic reconnection, buoyancy, flaring, wind mass 
loss, etc.) provide an interpretatlonal framewor)< for other astrophysica 1 
systems in which such processes are thought to occur (e.g., Rosner 1982). 

These major developments in stellar x-ray astronomy follow on the heels 
of a similarly major readjustment in our understanding of solar surface 
activity, which resulted primarily from the wealth of data gathered by the 
Slcy lab ATM Instruments over a decade ago (see reviews by Withbroe and Noyes 
1977; Vaiana and Rosner 1978: and the monographs of the Skyl ab Workshops, 
edited by Zlrker [1978], Sturrock [1980], and Orrall [1981]). Since then, it 


has become clear that the entire range of solar coronal activity is governed 
by the interaction between the surface magnetic fields of the Sun and the 
ambient plasma; thus, magnetic fields play a central role in the energetics of 
not only solar flares (as has been long realized), but also of the far more 
quiescent "Quiet Sun". The physical processes which control local conditions 

— hydrodynamics, heat conduction, radiation, mechanical heating by waves and 
currents — Involve the complications of microscopic plasma effects well-known 
to the laboratory plasma physicist, albeit under a vastly different 
combination of physical conditions and spatial scales. In addition, processes 
less familiar to the laboratory physicist -- stochastic magnetic dynamo 
activity, magnetic flux concentration by turbulent fluids, magnetic buoyancy 

— must be considered. These microscopic and macroscopic processes 
undoubtedly also play a central role in the physics of activity on other 
stars, as well as in the physics of rather different astrophysical systems, 
such as the interstellar medium, supernova remnants, accretion disks, and 
similar (but possibly still more "exotic") objects. In many of these cases, 
the lack of spatial resolution mitigates against direct confrontation of 
theories with observations; the astrophysicist must instead rely upon 
extrapolation from (analogous) cases in which theory and observations can be 
compared relatively directly. Thus, solar and stellar physics can play a 
central role in the interpretation of plasma phenomena in distant astronomical 
objects . 

The emergence of stellar x-ray astronomy us a new astrophysical 
discipline hai also played a central role In the recent resurgence of interest 
in the "solar-stellar connection". This resulted from the serendipitous 
conjunction of stellar x-ray observations provided by the now-defunct Ei nstei n 
Observatory (Vaiana et al . 1981; Helfand and Caillault 1982; Linsky 1982; 
Valana 1982, 1983; Golub 1983a, b; Haisch 1983; Johnson 1983a; Stern 1983; 
Serlo 1984), the renewed exploitation of ground-based Ca II observations 
(Wilson 1978; Vaughan et nl . 1982), the spectroscopic stellar observations 
from the International Ultraviolet Explorer (Dupree 1982, 1983; Linsky 1982, 
1983), and the increased sensitivity in radio observations made possible by 
the VLA (viz., Gibson 1981; Dulk 1985). This confluence has led to 
substantial synergism, in that methods of analysis first developed in one of 
these contexts have been rapidly exploited in the others; the studies of the 
variation of activity levels with various stellar attributes (such as 
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rotation) are the most prominent example of this synergism. 

In this review, we shall focus on the main physics problems which arise 
when confronting the stellar x-ray data. We shall not examine the details of 
the data, except insofar as such studies lead us directly to an important 
physical issue; we will Instead refer the reader to the specific literature 
(Including the large number of available reviews of the data) where such 
details are provided. A central issue throughout our review is the exte.it to 
which the extrapolation from solar physics is Justified, and to define (if 
possible) the limits to such extrapolation. This issue is of particular 
concern as one confronts x-ray emission from stars which are not traditionally 
thought to be "solar-like” (such as OB stars and T Tauri stars): while the 
stellar astronomer "bred" in the solar context may be tempted to apply the 
"solar analogy", those who have focussed on the physics of such stars may have 
radically-dif f erent notions as to the origin of the x-ray emission (even if 
its presence has come as a great surprise to such workers!). 


1 1 . THE BASIC PROPERTIES OF X- !^Y I SS ION FROM ST ARS 


Observations of quiescent x-ray emission from stars were first made in 
the middle 1970' s (Catura, Acton, and Johnson 1975; Mewe e t a l . 1975), and the 
EXOSAT satellite is currently acquiring both iiftaging and spectroscopic data on 
such emission (viz., Landinl et al . 1984a, b; Mewe 1984; Westergaard ez_^l_. 
1984). Nevei theless . virtually all of the published data to date on stellar 
x-ray emission ha^’e been obtained from two satellites, the first High Energy 

Ast.'onoray Observatory (H EAP 1) and the Einstein Observatory ( HEAP 

Arguably the mort dramatic impact has been made by images acquired from the 
latter instrument, and much of this review will focus on the observational and 
theoretical implications for our understanding of solar and stellar activity 
which have emerged from these results. Because of the wealth of reviews in 
the literature which summarize Uie observational details (Rosner and Vaiana 
1979; Linsky 1980. 1982, 1983a, b; Vaiana 1980, 1981a, b. 1982; Mewe. 
Schrijver, and Zwaan 1981; Golub 1983a. b; Johnson 1983a; Stern 1983; 
Cassinelli and McGregor 1984; Serio 1984), we will not emphasize the minutiae 
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of the observations, but rather focus on the major observational questions, 
the outstanding observational problems, and the theoretical implications of 
the data now In hand. 

What ought to be expected from stellar x-ray observations If the Sun is a 
prototypical late-type dwarf star? The visually dramatic charge In the 
observed solar activity level seen during the declining phase of the solar 
activity cycle by the ^y_lai>. so. t x-ray telescopes corresponds to a change of 
well below an order of magnitude in energy output at soft C 1/4 keV) x-ray 
wavelengths; and if one takes the full variation of coronal radiative output 
(excluding flares) during the course of an entire activity cycle into account, 
then one would expect roughly a factor of ten variation in the stellar soft 
x-ray luminosity (Pallavicinl e^jl.. 1981). As is now well-known, this is not 
at all what is observed; if the Sun were really prototypical, the E instein 
Observatory stellar surveys would have been meager affairs. Instead, the Sun 
lies near the bottom end of the observed range of x-ray luminosities of 
late-type stars, which spans over three decades. This fundamental result is 
vividly displayed by the stellar x-ray luminosity functions for late spectral 
type dwarf stars discussed below: clearly, the Sup's coronal luminosity is 

very modest when compared to that of other stars. 

The basic observational facts regarding the level of x-r.iy emission are 
simply summarized in the H-K diagram of Figure 1, and by the cumulative x-ray 
luminosity functions displayed in Figure 2 (for the basic observational data, 
see references to the surveys and observational reviews cited above, as well 
as those cited specifically below): 

(1) The data are consistent with aXL^^srf stars of spectral type dF 

through dM beit.g x-ray emitters, with quiescent luminosities ranging between 
26 31 -1 

roughly 10 and 10 ergs sec 

(2) stars earlier than (roughly) B5 are x-ray emitters, at emission 

29 34 -1 

levels ranging between 10 to 10 ergs sec , with a roughly constant 

-7 

x-ray-to-bolometrlc luminosity ratio of 10 (Hamden et_al^. 1979; Seward et 
al . 1979; Long and White 1980), independent of luminosity class. 


(3) There is a narrow spectral range — from roughly B8 to A5 -- on the 
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■ain sequence for which there is no credible evidence for any x-ray emission 
from spectroscopically normal stars and Am stars (Golub et al . 1982; Schmitt 
et a i . 1984): however, there is evidence for x-ray emission from Ap stars (see 
also Cash and Snow 1982). 

(4) Late spectral type giants and supergiants experience a cutoff in 
x-ray emission levels as one moves to later spectral type (Ayres and Linsky 
1980; Ayres e t al . 1981: Vaiana et_al_. 1981; Haisch and Simon 1982): the 
nature and location of this cutoff is discussed further below. 

(5) Pre-main sequence stars are very vigorous x-ray sources, and show a 
definite trend of decreasing mean emission levels with nominal stellar age 
(Felgelson 1983; Mag:rlo et a l . 1984: Nicela. Sciortino. and Serio 1984; Micela 
et al . 1985: Caillault and Helfand 1985). 

(6) The x-ray spectia of late-type stars are therral. with 

6 7 

single-component temperatures ranging between 10 to several times 10 K; 
however, there is considerable evidence that a single-component analysis is 
not an adequate description of spectra for sources with sufficient count 
statistics (Holt et al . 1979; Swank et al . 1981; Mewe ^_al_. 1982: Swank and 
Johnson 1982: Vedder ard Canizares 1983; Schrljver, Mewe, and Walter 1984), 
and that a continuous emission measure distribution in temperature is far mure 
appropriate (Majer et a l . 1984; Mewe '984: Schmitt 1984a; Schmitt et al . 
1984a) . 

(7) As would be expected from soiar-lixe activity behavior, late-spectral 
type main sequence stars (Haisch et a i . 1980, 1983: Kahler et_al_. 1982) and 
evolved stars in close binary systems (e.g.. the RS CVn's; viz ., Walter et a l . 
1980; Agrawal et a 1 . 1983) show episodic highly transient behavior in x-rays: 
Indeed, very young main sequence and nre-main sequence stars show the most 
dramatic range of such behavior (Stern 1982: Montmerle e t a 1 . 1983; see 
especially Stern, Underwood, and Antiochos 1983 for com[iar isens ) . In 
contrast, early spectral type stars show little evidence for such short time 
scale variability, but do show evidence for longer-period fluctuations (e.g.. 
Coliura et a i . 1984). We shall not consider variability, especially of the 
flare or transient type, any further in this review, and Instead refer the 
reader to the reviews of Haisch (1983), Gibson (1983), Charles (1983), and 
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Mullen (1983). 

(8) X-ray emission levels from close binary systems, such hs the RS CVn 
(viz., Walter et a l . 1980) and the W UMa stars (viz.. Crudd-ice and Dupree 
1984), in which accretion is not thou(;ht to play a major role, are 
nevertheless substantially enhanced over that from effecLivniy single stars 
whose classical stellar characteristics (e.g,. mass, radius, effective 
temperature) are similar to those of the individual binary components (see 
reviews by Charles 1983: Dupree 1983; Linsl<y 1984). 


III. THF. SOLAR ANALOGY AND THE M ODELING OF X-RAY EM 1 SSI O N FR OM LATE- TYP E 
STARS 

The first issue to be dealt with is how the x-ray data just summarized 
are to be interpreted. If we confine our attention to the observations of 
late spectral type stars, then the obvious candidate analogue * *he soiar 
corona: however, if such extrapolation is to make sense at all, • fs it ought 
to be the case that canonical stellar activity indicators other ^ lian x-ray 
emission have plausible counterparts in the solar data; and that, conversely, 
the kinds of correlations between such indicators seen on the Sun are also 
observed to hold in the stellar data. 

Perhaps the best place to start is with d brief remlndei of what the Sun 
looks like when viewed at a variety of wavelengths; the extensive full-disk 
x-ray image catalogue of Zombeck e^^ . (1979) gives perhaps the best overview 
to date. Prom such images, it Is evident that -- even on the immediate, 
qualitative level — the Intensity of x-ray emission is very inhomogeneous, 
and correlates rather well spatially with emission from lower-temperature 
species, as well as with the photospheric magnetic field Is^alana aid Rosner 
1978 and references therein). Furthermore, it is evident that the spatial 
contrast of activity levels becomes greater as the observational wavelength 
decreases: minimal in the optical, and maximal in the x-ray (the extremes -- 
hard x-rays and gamma rays — are extremely localized, and are seen virtually 
only during the course of solar flares). This variation of contrast with 
wavelength is Illustrated by the composite shown ns Fig. 6 In Valana and 
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Rosner (1978), and allows us to conclude that the activity contrast is 
greatest at precisely those wavelengths at which most of the radiative 
emission from a hot plasma — the chromosphere and corona — occurs. 
Unfortunately, the terrestrial atmosphere is strongly absorbing at these 
wavelengths, and it is therefore not totally surprising that, prior to about 
1975, there were n^ direct observations of steady coronal activity from 
solar-like stars (e.g., Mewe 1979; Ulmschneider 1979): it was not unt'l 
roughly after that date that space-borne telescopes of sufficient sensitivity 
at short wavelengths became available. Until then, studies of stellar 


1 

This difficulty was of course readily circumvented (at least partially) by 
finding signatures of solar surface activity which (i) can be seen at 
wavelengths accessible to ground-based observers, and (ii) can be detected in 
spatially-integrated light (i.e., if one views the Sun as a stellar-like point 
source; viz., Skumanich e t_a 1 . 1984). If such signatures exist, then, as 
pointed out above, it is in principle possible to test whether stars show the 
kinds of correlations between activity Indicators observed on the Sun: 
sensible extrapolation from solar physics hinges on a positive result. 
Indeed, a variety of such indicators exist, ranging from modulation of the 
integrated stellar light in the visible (due to "spots"; e.g., Oskanyan et al . 
1977; Willson et al . 1981) to emission in the Ca II H and K lines (Wilson 
1966, 1978; Vaughan 1983) and in He 10830 A (Zirin 1975, 1982). The observed 
correlations between these ground-based activity indicators are what one would 
expect on the basis of solar observations, and so the assumption that what we 
are observing indeed can be understood by extrapolation from the Sun appears 
to be on solid ground. More recently, these ground-based studies have 
received renewed impetus because of the exploitation of high-resolution 
optical spectroscopy to detect the Zeeman splitting of magnetical ly-sensitive 
lines (Robinson, Worden, and Harvey 1980; Marcy 1983), allowing the kind of 
correlative studies of coronal and photospheric activity which have been 
successfully carried out for the Sun (viz., Vaughan 1984). 


activity relied upon ground-based proxy indicators for stellar activity — 
modulation of integrated light, Ca 1 1 H and K emission, and observations of 
the He 10830 line — which indeed strongly supported the very attractive 
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proposition that solar-like surface activity is not peculiar to the Sun alone. 
However, It must be kept in mind that ground-based diagnostic tools share the 
common defect that their very visibility from ground is evidence for rather 
non-solar behavior: that is. if we were to view the Sun as a star, and place 

it at typical stellar distances, little regarding the properties of solar 
activity could be deduced on the basis of these three classical ground-based 
diagnostics — solar activity is simply not sufficiently vigorous. The use of 
ground-based diagnostics has thus the effect of selecting a very special 
sub-class of stars, namely those which show presumably solar- like activity 
behaviour, but to an extreme: and we therefore run the substantial risk of 

biasing our conclusions regarding the general behavior of stellar surface 
activity . 

These kinds of selection effects ought to be avoidable if the sample of 

stars considered is volume-limited, rather than flux-limited. In order to 

maximize the size of the sample, one ought therefore observe stellar activity 

at those wavelengths at which the hot plasma emits most of its energy: and, 

with sufficient sensitivity, one is then in a position to explore a 

significant portion of the coronal iumino«’ty function (i.e., the probability 

2 

distribution function tor the coronal luminosity) . In this way, one might 


2 

A crucial aside is that in conducting such volume-limited surveys, it is 
essential to consider both the detections and the upper limits (e.g., failures 
to detect); one can then use so-called "detections-and-bounds" (Avni e t al . 
1980: see Rosner 1983 for details on applications to stellar observations) or 
"survival analysis" (Schmitt, Golub, Hamden, and Rosner 1984: Schmitt 1984b; 
Feigelson and Nelson 1985; Schmitt and Rosner 1985) methods in order to 
construct maximum likelihood distribution functions for stellar x-ray 
parameters. In particular, Schmitt, Golub, Hamden, and Rosner (1984) have 
shown by example that disregard of the detection upper limits can lead to 
conclusions entirely opposite to those reached by considering the full set of 
information available. 


hope to observe the full range of possible stellar activity levels, and 
minimize selection effects. Thus, if the Sun is to be any guide, one ought to 
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observe stellar activity most optimally at UV and shorter wavelengths. In 
saying this, one must of course be prepared for the possibility that plasma 
temperature regimes associated with stellar activity are ratlier different than 
is the case for the Sun: the difference in solar coronal temperatures of the 
Quiet Sun and coronal holes as deduced by, for example, soft x-ray diagnostics 
(e.g. , Maxson and Vaiana 1977) should stand as an important warning. Indeed, 
definition of the domain of applicability of the "magnetically-confined" 
corona prototype is itself a central problem of stellar activity modeling, 
particularly for late-type giants and supergiants (viz., Dupree 1982; Linsky 
1982) . 

Let us then assume that we have obtained a relatively bias-free set of 
observations of stellar coronal emission, and ask how the data are to be 
interpreted. If we again confine our attention to stars of 1 ite spectral 
type, then the natural starting point is the basic fact of the structuring of 
the (solar) corona. The existence of coronal structures, and the important 
role played by magnetic fields in such structuring (which is most obvious in 
solar active regions) was recognized long before high resolution EUV and x-ray 
images were available (see, for example, Billings and Saito 1964). However, 
the availability of on-disK observations of the corona in the EUV and in 
x-rays revealed the ubiquitous loop structuring of the emitting regions 
throughout the solar corona, as well as the close connection with regions of 
strong surface magnetic fields (Vaiana, Krieger, and Van Speybroeck 1970). A 
series of rocket flights led to semi -empirical modelling of coronal loop 
structures and coronal magnetic fields (Vaiana, Krieger. and Timothy 1973); 
and the long duration Skylab missions permitted study of the temporal 
evolution of the emission regions in the corona, and of the quantitative link 
between magnetic fields and coronal emission. 

These data show the solar corona to be the superposition of a variety of 
coronal structures, whose physical parameters (size, temperature, density) 
cover a wide range (see Withbroe and Noyes 1977), and whose relation to the 
topology of the defining magnetic fields is quite strict; virtually all of the 
coronal emission comes from topologically closed magnetic structures (Van 
Speybroeck, Krieger, and Vaiana 1970; Vaiana, Krieger, and Timothy 1973; 
McIntosh et al_. 1976), with regions of topologically open magnetic fields 


characterized by severely lowered x-ray fluxes (Maxson and Vaiana 1977). In 


all of the following discussion, we will assume that if the solar analogy 
applies, then the high temperature plasma seen in x-ray emission from stars is 
in fact predominantly confined in such topologically closed loop structures. 

The importance of structuring in coronal physics was put on a firm 
theoretical foundation in the late 1970's (Rosner, Tucker and Vaiana 1978: 
Rosner and Vaiana 1978; Craig, McClymont, and Underwood 1978; Vesecky et al . 
1979; see also Priest 1981, Withbroe 1981, and Hammer 1983); the key 
realizations were, first, that the coronal magnetic fields associated with 
individual loop structures physically confined the plasma (so that mass loss 
to the solar wind is insignificant for these structures) and, second, that 
these loops could be considered as relatively independent "mini-coronae" 
because of the thermally-insulating effect of the confining magnetic fields 
(Landini and Monsignor i-Fossi 1975). Thus, the solar corona can be viewed — 
at any given time — as an ensemble of these building blocks whose precise 
statistical mix of physical properties (size scale, mean field strength, mean 
plasma temperature and pressure) is determined by the stochastic eruption of 
magnetic flux from the solar interior (Vaiana and Rosner 1978). That is. the 
emerging solar surface fields most certainly play a passive role in 
determining the coronal topology (and thus determine the confinement 
properties of the coronal gas): in addition, they may also play an active role 
by providing the direct means for dissipation of mechanical energy flux to 
heat the coronal plasma in_ (Gold 1964; Tucker 1973; Rosner. Tucker, and 
Vaiana 1978; lonson 1978; see reviews by Wentzel 1978 and Kuperus, lonson, and 
Spicer 1981, and the more recent work of Parker 1983a, b). These various 
considerations have a remarkably simple consequence to lowest order in the 
physical description of such structures: if we consider the balance between 
only local heating (by whatever means), thermal conduction, and (optically- 
thin) radiatiative losses in a hydrostatic (confined) plasma, one directly 
obtains a testa ble scaling relation involving the observable quantities loop 
length, (peak) coronal temperature, and loop pressure (Rosner. Tucker, and 
Vaiana 1978; Craig, McClymont, and Underwood 1978). In concert with more 
sophisticated observations, these scaling relations can be adumbrated to 
include more physics, such as gravitational stratification (Serio et al . 1981; 
Wragg and Priest 1981) and steady flows (viz., review by Priest 1981), but the 
essential element remains: the coronal plasma is a bright source of x-ray 
emission because of magnetic confinement. Thus, the first point to recognize 
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when discussing observations of x-ray emission from solar-like stars is that 
the observed emission very likely derives from a statistical mix of 
structures, so that any attempts to describe the observations by appealing to 
a "mean" atmosphere must be viewed with some trepidation. 

Given the extremely complicated mixture of loop structures which are 

found in the solar corona, can spatially unresolved stellar observations 

(which necessarily average over the entire ensemble of structures in a stellar 

atmosphere) be of any real use in studying the nature of stellar activity and. 

more specifically, the coronal heating mechanism? Here nature has been kind 

to us, in that it is now clear that variability in emission properties is a 

common feature of stellar x-ray emission; the hope is then to use the 

variability as a probe of the kind of structuring to be expected in stellar 
3 

coronae. The observed variability is of course not surprising, since the 


3 

The outstanding exception to this is the use of eclipses in relatively close 
binary systems to disentangle the geometric relationship between different 
coronal components, as has been illustrated by the study of the eclipsing KS 
CVn system AR Lac by Walter, Gibson, and Basri (1983; see also White et al. 
1984) . 


solar corona is known to be variable on a remarkable wide range of time scales 
(even when flares are excluded), as is wbll-illustrated in the sequential 
x-ray images of the Sun published by Zombeck et al . (1979): these sequences 
show five solar rotations observed from SkyJaJ). and demonstrate that on 
essentially any time scale of interest — from minutes, hours, and days to 
months -- the Sun may be regarded as a variable x-ray source. These 
variations have a wide dynamic range: thus, an example of a rather dramatic 
variation Is seen during the fourth and fifth rotations, when the corona 
varied from a solar maximum-like configuration to one more characteristic of a 
solar minimum-like configuration during a time span of less than seven days. 
However the full variation in coronal composition during different portions 
of the solar cycle is even more extreme than that: for example, x-ray data 
obtained near solar minimum show an entirely different kind of corona, winch 
is dominated by large numbers of small active regions (the so-called 
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"ephemeral active regions"). These snail regions are short-lived, do not 
develop sunspots and have relatively low x-ray emission measure, but are found 
to contain as much magnetic flux as the larger, longer-lived active regions 
found at times of higher solar activity (Golub 1980). 

These solar coronal variations can be quantified by representing the 

variation of the total emission by a superposition of the various known 

atmospheric components, and by accounting both for the solar cycle variation 

of each component, and for the evolutionary histories of the different types 

of emerging field regions associated with these coronal plasma constituents. 

From quantitative studies of this kind, Golub (1983) concluded that the 

average solar corona is dominated by a single type of structure approximately 

half of the time. That is, at solar maximum, the entire coronal x-ray 

emission can be characterized by the properties of large active regions, with 

0 0 

a temperature (T) of approximately 10 ' k and an emission measure (EM) of 
, 50.5 -3 

approximately 10 cm At solar minimum, the corona is characterized by 

the small emission regions (together with their diffused — evolved — 

6.3 49 -3 

states), with T approximately 10 K and EM approximately 10 cm This 

then suggests that It may be possible to describe the morphology of stellar 
coronal emission by superposition of a minimal number of distinct emission 
components, especially if models are further constrained by combining the 
x-ray data with contemporaneous observations of the cooler underlying 
atmosphere (cf., Mewe and Zwaan 1980; Ayres, Marstad, and Linsky 1981; 

Pallavicini et al . 1982; Schrijver et al . ^983; Giampapa et al . 1985). As we 

shall discuss below, there is in fact some evidence that such a minimal 

multi-component description is an apt summary of the stellar x-ray data. 

These preliminaries allow us to place x-ray observations of late-type 
stars into an interpretational framework which is well-tested. However, it 
may well be that this (solar) analogy is Inappropriate, and hence it is 
essentia] to continually devise observational tests to verify its 
applicaoility . The most prominent example of how the exact solar analogy may 
fall for a "solar-type" star is AR Lac: Walter, Gibson, and Basri (1983) have 

in fact shown that significant x-ray emission comes from regions well-removed 
from the stellar surfaces of the two binary companions ii; *'iis system, so that 
a model of "interacting loops" (In which the magnetospheres of the two 
companions Interact to produce local heating, and hence x-ray emission; Uchida 
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and Sakurai 1983; Uchida 1984) may be a far more appropriate description of 
the source morphology. 


I V . MODELING X-RA Y EMI SSIO N FROM EAR LY-TYP E STAR S 

In contrast to the case of late-type stars, we do not have ready access 
to a nearby astronomical object which shares its stellar properties with that 
of the hot stars; it is therefore not entirely surprising that the relative 
state of our understanding of the outer atmospheres of early spectral type 
stars is quite primitive when compared with that of solar-type stars. Since 
this subject area is hardly mature, we will confine ourselves to a very brief 
review of the principal physics issues: and refer the reader to the recent 
reviews by Cassinelli (1984) and Hearn (1984) for further details. 

Studies of possible x-ray emission from the atmospheres of early-type 
stars proper (as opposed to, for example, the terminal shocks associated with 
their winds) which proceeded the definitive detection of soft x-rays from 
these stars (Hamden ot aj_. 1979: Seward e^aj. 1979) focussed on models which 
were essentially analogues of the solar atmosphere in extremis: a hot corona 
near the stellar surface, underlying a far cooler high-speed wind (Cassinelli 
and 0*.sen 1979). It is by now apparent that this kind of model, as well as 
its various variani's (Rosner and Vaiana 1980; Stewart and Fabian 1981; 
Cassinelli 1984: Uchida 1984: Waldron 1984; Underhill 1984) cannot easily 
reconcile the observed extinction (at energies below roughly 0.5 keV) and the 
spectrum (whose derivation is dominated by photons with energies above 
approximately 0.5 keV) derived from the x-ray data with the imputed wind mass 
column densities derived from radio and IR observations (see Cassinelli e t a 1 . 
1981 for details). That is not to say that, by Invoking sufficient complexity 
for this class of models, an acceptable fit to the data cannot be achieved; 
instead, the real difficulty is to account in a natural and plausible way for 
the presence of any plasma heating in th e first place -- to appeal to, for 
example, the solar paradigm of magnetic and/or wave heating is not easily 
physically motivated, and hardly accounts in a natural way for one of the 
principal observational constraints. namely the remarkably constant 


x-ray-to-bolometric luminosity ratio as a function of spectral type. 
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An alternative physical picture of what might be occuring was proposed by 
Lucy and White (1980; also Lucy 1982), who revived earlier suggestions by Lucy 
and Solomon (1970), Hearn (1975), Nelson and Hearn (1978), and McGregor, 
Hartmann, and Raymond (1979) that the massive radiatively-driven wind might be 
unstable (see also Carlberg 1980; Kahn 1981; Owocki and Rybicki 1984), and 
argued that the higher-density "blobs" formed as a result of the instability 
would lead to the formation of shocks (and associated high-energy radiative 
emission) throughout the unstable region of the wind outflow. The advantages 
of this model are several; (1) it is now Indisputable that the 
radiatively-driven winds of OB stars are Indeed linearly unstable to density 
perturbations (thus overcoming a principal objection to the corona-wind 
model); (ii) shock-heated matter is expected to be distributed throughout the 
wind, thus substantially easing the difficulty of reconciling the observed 
wind mass column aensities and soft x-ray extinction; (iii) the emission is 
expected to be naturally time-variable, as it is the superposition of emission 
from many distinct shocked regions of the wind. This latter expectation is 
indeed in accord with observations of variability at x-ray wavelengths 
(Collura et al . 1984). 

However, we do not believe that this subject is as yet closed because 
several major gaps still separate the observations from a reasonably-complete 
model : 

(i) Although linear Instability has been' demonstrated, the non-linear 
development of the instability to the point of "blobs" (or their effective 
equivalent) remains to be studied; 

(ii) The simple phenomenological shock model (in which the shock strength 
is fixed) does not in fact give a good account of the x-ray data (Cassinelli 
and Swank 1983); what is required is a distribution of shock strengths, which 
must at this juncture remain entirely ad hoc (because no theory of shock 
formation, which evolves the flow non-linearly from the linear domain, exists 
at present ) . 

Resolution of these various issues will unquestionably involve (in 
addition to the theoretical work just alluded to) a substantial effort at 


-16- 


mult l-wavelength observations (see, for example, Caillault et al . 1985 and 
White and Becker 1983), preferably carried out with (near-) simultaneity in 
order to correlate emission levels and variability over the entire range of 
wavelengths in which non-thermal emission has been observed. In addition, the 
Einstein observations provided a tantalizing hint that higher spectral 
resolution and sensitivity (especially below 1 keV) could provide decisive 
tests of extant models (viz., Cassinelll and Swank 1983). 


V. THE PHYSICS OF STELLAR X-RAY EMISSION 


The most obvious question which the data discussed above engenders is 
simply why stars emit x-rays. The answer to this question is not 
straightforward, and indeed for the early spectral type stars cannot as yet be 
reliably answered. Our understanding of the x-ray emission processes for the 
late spectral type stars is heavily dependent upon the applicability of the 
solar analogy discussed in the previous section; and, in so far as this 
analogy applies, we do understand much of the basic physics leading to stellar 
x-ray emission. 

( i ) Corre lati on Betwee n Stel lar X-ray Emiss ion and Other Stellar 

Attributes . 


The first question is simply observational: given the large dispersion in 
emission levels for stars of any given spectral type and luminosity class 
(Figure 2), can we determine which physical attributes of stars account for 
the observed variance? Thus, because the total x-ray luminosity appears to be 
only very weakly related to the effective surface temperature (and hence to 
the level of surface fluid turbulence) for late-type main sequence stars, some 
stellar properties which do not figure in placing a star in the H-R diagi-am 
must determine the actual emission levels (this is not the case for the early 
spectral type stars, as discussed below). 

As posed, this issue is a classical statistics problem (common factor 

analysis), with the not inconsiderable complication that the data is heavily 

4 

censored (e.g., contains many upper bounds). What these missing parameters 
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Indeed, some work in applying classical common factor analysis to data which 
are only weakly censored (and so for which the upper limits — which were 
ignored — carry only a very limited amount of information) has been carried 
out by C. J. Schrijver (1982), who showed that for the limited sample he 
considered, stellar rotation was the dominant factor in determining the level 
of stellar x-ray emission; see also Majer et al . (1984) in the context of the 
"rotation-activity connection" for RS CVn stars. 


might be was suggested early on by HEAP 1 observations (Walter et al . 1980) 
and subsequent E instein Observatory surveys (Walter and Boywer 1981) of RS CVn 
stars, which suggested that rotation is a significant determinant of coronal 
luminosity for these binary systems: we shall return to this particular issue 
below. Such a connection between rotation and chromospheric activity level 
was long known from Ca II data (Kraft 19S7); and discovery of a coronal 
luminosity-rotation correlation for Isolated and effectively single late-type 
dwarf stars followed from analysis of the Einstei n/Center for Astrophysics 
(CfA) stellar survey data (Pallavicini et al . 1981, 1982). In contrast, it 
was evident from the first observations of x-ray emission from early spectral 
type stars that the dominant determinant of emission levels was the bolometrlc 
luminosity (Hamden et al . 1979; Seward et al . 1979; Long and White 1980); 
thus, for these stars, the location in the H-R diagram is Indeed a very good 
predictor of x-ray emission levels, 

( a ) Stellar x-ray emission and stellar rotation. 

The basic points regarding the correlation of x-ray emission levels with 
rotation rate are summarized in Figure 4, which represents an update of the 
data presented by Pallavicini et al . (1981): the rotation dependence of 
coronal emission for stars of spectral type F7 to M5, and luminosity classes 
III, IV and V, is quite clear. The linea.' bes.' !• to the total data sample 
is based on detections alone, and indicates ..ut L is proportional to Die 

X 

square of the stellar rotation rate. In contrast, the early-type stars do not 
show any rotation dependence of x-ruy emission levels (Pallavicini et al . 
1981), Some controversy has occurred regarding the precise slope of the power 



- 18 - 


law relation connecting ^he rotation period (cf., Walter and Boywer 

1981; Pallavlclnl et al ■ 1981, 1982; Rengarajan and Verma 1983; Rengarajan 

1984), but it appears that the fundamental difficulty (which remains in place 
to date) is the introduction of possibly large biases because of the ill- 
defined nature of the various samples of stars for which both rotation rates 
and x-ray luminosities are known. We also note that similar difficulties 
beset determination of a correlation between the stellar activity level and 
characteristic Rossby number (viz., Noyes et al . 1984), discussed further in 

the next section: as emphasized by Schmitt, Golub, Hamden, and Rosner (1985), 
such correlations using x-ray data must be very cautiously interpreted, 
especially if no account is taken of the presence of upper bounds in 
measurement of x-ray luminosity, stellar rotation rate, or both (cf., Mangenay 
and Praderie 1984; Marilli and Catalano 1984). 

For one class of late-type stars, it has however been possible to conduct 

a systematic study of the relation of stellar rotation to surface activity: 

the RS CVn stars. These stars form a class cf binaries whose components tend 

to be somewhat evolved; their orbital periods are generally less than about 14 

days (Hall 1976), although it is common to include much longer-period systems 

in this category (such as Capella). By comparing the x-ray to bolometrlc 

luminosity ratio L /l with the rotation periods of these systems (which 

X bo 1 

correspond to the orbital periods in the synchronous shorter-period systems), 

- 1.2 

Walter and Boywer (1980) showed that this ratio varied as (period) . The 
class of stars examined is however quite heterogeneous, and one must therefore 
ask whether some other stellar attribute could possibly lead to this 
correlation: this possibility was raised by Pallavicini et al . (1981), who 
pointed out that L was not correlated with period for this sample of stars, 

X 

but that the bolometrlc luminosity w^. This problem has been studied in 

detail by Rengarajan and Verma (1983), who have shown that the bolometrlc 

0.7 

luminosity is Indeed related to the orbital period as (period) for a much 

larger sample of F, G, and K eclipsing binaries (cf.. Young and Konlges 1977), 
1 .24 

and as (period) for the more restricted Waiter and Boywer sample: hence it 
appears that essentially all of the correlation between L /l and period is 

X bol 

due to the correlation between L ^^d period, a conclusion which has been 

bol 

further strengthened by a formal principal component analysis of the RS CVn 
data which explicitly shows that there is no significant "rotation-activity" 
correlation for these stars (Majer et a l . 1984). Indeed, the correlation 
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between and period is to be expected bcccjse of the simple fact that for 
tldally-coupled close binaries. the stellar radius Is roughly proportional to 
the stellar separation, so that — via Kepler's law — there Is a correlation 
between bolometrlc luminosity and orbital (and hence rotation) period; see 
Young and Konlges (1977). 


(b) Stellar x-ra y emission a.hd-Si®Uar age • 

A most basic question Is whether stellar activity levels show any 
evolutionary effects as a given star moves along its track in the H-R diagram, 
from protostellar object of well-evolved giant or supergiant. This question 
has been systematically addressed by several extensive surveys of stellar 
clusters and associations, ranging from T Taurl associations (Ku and Chanan 
1979; Peigelson and Kriss 1981; Montmerle e^jal,. 1983) to the classic nearby 
young open clusters Hyades (Stern et al . 1981; Zolcinskl et al . 1982), 

Pleiades (Helfand and Caillault 1982; Hlcela, Sciortlno. and Serio 1984; 
Mlcela et a l . 1985; Caillault and Helfand 1985), and Ursa Majoris (Walter e^ 

1984). 

Let us begin by considering the youngest stellar objects. During the 

final phases of star formation, a contracting protostar is enveloped in a 

shroud of accreting gas, so that until recently, this phase of stellar birth 

was invisible to us (viz., Gertz, Black, and Solomon 1984). However, It had 

been speculated by theorists that such objects might be x-ray sources, and 

hence observable in spite of the large obscuration in the optical: thus, 

Bisnovatyi-Kogan and Lamzin (1977) suggested that if the outer atmospheres of 

T Taurl stars were expanding outward (e.g., in the form of a stellar wind), 

and are thermally-driven, then a necessary consequence would be the presence 

of a hot (million degree) corona near the base of the outflow region; the 

34 - 1 

corresponding x-ray luminosity would then be of order 10 erg sec Further 
x-ray emission would be also expected from the collision between this outflow 
and the ambient interstellar matter (viz,, Schwartz 1978). In contrast, if 
the atmospheres of such stars were dominated by accretion from the ambient 
interstellar matter, tlien one might expect emission due to the impact of the 
accreting matter on the stellar surface (viz., Ulrich 1976; Mundt 1981), with 
the temperature of the impact-heated matter largely determined by the infall 
speed. The difficulty with these various predictions is that they were 
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largely contradicted by the earliest observational studies: thus, observations 
by ANS and SAS-3 resulted in the detection of only one source which co’ild be 
associated with one cr more pre-main sequence objects (den Boggen'^e et al . 
1978; Bradt and Kelley 1979), but whose luminosity was clearly far below that 
supposed by Bisnovatyi-Kognn and Lamzin (1977), and whose source spectrum was 
quite a bit harder than one would expect for the accretion models. It is from 
this historical perspective that the striking Einstei n Observatory data from 
the Orion Nebula (Kii and Chanan 1979: Ku et al . 1982), the Taurus-Aurlga 
complex (Gahm 1980; Feigelson and DeCampli 1981; Walter and Kuhl 1981, 1983), 
the Rho Ophiuchl cloud (Monlmerie et__^. 1983), and the Chamaeleon cloud 
(Feigelson and Kriss 1981) must be viewed. 

These Ei nste in observations showed, first, that pre-main sequence stars 
without question constitute a new class of soft x-ray sources: second, that 
their emission — when seen — is of order 1,000 times more vigorous than that 
of main sequence stars of similar spectral type; third, that their emission is 
extremely variable; fourth, that the level of x-ray emission is correlated 
with the stellar optical magnitude, but not with the classic indicators of 
pre-main sequence "activity" (such as H alpha emission): and. fifth, that the 
spectrum is cut-off at low energies (as is expected from absorption by the 
intervening mutter within the star-forming region), and seems to be 
intrinsically hard (for a stellar source). These observations, which are 
reviewed in much greater detail by Feigelson (1983), raise a number of 
immediate questions: 

I 

Wh at is the origin of th e x-ray emission? The basic question to be 

answered is how apt the "solar analogy" really is for these stars. Can one 
regard the observed x-rays as the radiative emission from solar-like activity 
centers, or is the observed emission more properly regarded as the entire 
consequence of continual "flaring"? That is. while there is no doubt that 
substantial transient activity does take place (most beautifully shown by the 
results of Montmerle et_al_. 1983), the»'e is considerable question whether a 
less-variable "background" atmosphere Is present. It is especially intiiguing 
that this notion of continual very small-scale flaring as lh^e_ ul t imate coronal 
heating process has been again raised as a possibility for the Sun itself 
(viz., bin et_al_. 1984). In order to answer these questions in the case of 
very young stars, it is necessary to (i) place strong constraints on the x-ray 
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spectrum of the observed emission (In order to isolate the relative locations 
of the x-ray emitting matter and the cooler plasma associated with prc-mrln 
sequence stellar surface activity); (11) determine the extent of va''<eblllty 
(in order to determine whether flaring can entirely account for the obser''ed 
emission); and (ill) determine a bound on the level of "quiescent" emls».ion 
which Is comparable to that of more evolved late-type stars. 

What is the ulti mate cause of the activity ? Closely related to the above 
question is the problem of accounting for the underlying cause of "activity" 
stars. If. as in the case of late-type stars, the activity (whether transient 
or not) is the result of magnetic activity on the surface of these stars, then 
one would expert a similar sensitivity to the characteristics of the rotation 
characteristics and the convective interiors of these stars. In order to test 
these expectations, it is necessary to establish the x-ray luminosity 
functions for these very young starsi in a variety of star formation regions, 
••'’d then to explore their correlations with the rotation rale distribution, 
spectral type distributions, and age of the cluster stars; this program 
remains for the future, so that this question must remain unanswered at 
present . 


What happens once stars pass out of the protostel lar phase? A most 
useful comparison involves observations of the two classic nearby open 
clusters, the Pleiades (Caillault and Helfand 1985; Micela et al . 1984, 1985) 
and the Hyades (Stern et a 1 . 1981; Westergaard et a 1 . 1984) by the E inste in 
Observatory. Figure 4 gives the basic result: comparison of the cumulative 
x-ray luminosity functions for G stars in the Pleiades and Hyades open 
clusters, and for field G stars, shows a clear decline in activity levels with 
age. This result is as one would expect on the basis of the by now classical 
picture of the evolution of stellar activity: as the star evolves, its surface 
activity leads to spindown, and hence reduced surface activity, lowered 
spindown rate, and so forth (iUermann 1946; Schwarzschl Id 1948; Schat/man 
1962; S)(umanich 1972; Kraft 1967). However, this scheme is complicated by the 
fact that the data argue against a simple monotonic decline in activity 
levels. The particular problem is raised by the relative inactivity of both 
the K and M stars in the Pleiades; Indeed, there are no definitive detections 
of M dwarf stars in this open cluster. This puzzle is compounded by the fact 
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that the rotation rates of the Pleiades K stars are anoaalously high (Stauffer 
et al . 1984). 

Stauffer e t al . and Micela et_al_. (1985) argue that this anomaly can be 
understood by tracing the detailed evolution of the Internal rotation profile 
of these stars as they descend to the main sequence. That is, rapid evolution 
on the radiative track should lead to spin-up of the e ntire star (with modest 
radial differential rotation), followed by rapid spindown (by magnetic 
braking) of the outer (visible) convecting envelope (leading to significant 
differential rotation between the outer envelope and the radiative core; cf , 
Endal and Sofia 1981, and especially Tassoul and Tassoul 1984); the internal 
rotation profile in the latter case is particularly suited for vigorous 
magnetic dynamo activity (Rosner 1980; Spiegel and Weiss 1980; Golub e t al . 
1981; Schussler 1983). Hence, since the Pleiades K stars are thought to fall 
in the former category, and the G stars in this cluster in the latter 
category, the dominance of G star activity can be understood. 

( c ) Stellar x-ray emi ssion a nd sp e ctral type . 

For stars ranging from spectral type early 0 to roughly B5, the level of 

x-ray emission is proportional to the stars' bolometric luminosity; the 

-7 

proportionality constant is approximately 10 (see Pallavicini et aj_. 1981, 

and references therein). A transition between this "early spectral type" 
behavior, and the behavior of later spectral type stars (whose emission levels 

dominantly depend on rotation rate) occurs in the mid-A spectral type range 

(Schmitt et al . 1984); the coincidence between the location of this transition 
and the location of the onset of surface convection leads to the obvious 
suggestion that surface convection is indeed the necessary prerequisite for 
solar-like stellar activity (see also Bohm-Vitense and Dettmann 1980). 

Other than the aforementioned transition stage from "early-type" 
behavior, late spectral type stars generally do not show any evidence for a 
primary dependence of x-ray emission levels on spectral type; this fact is 
best illustrated by the sequence of x-ray luminosity functions for distinct 
spectral types shown in Figure 2. However, a major exception to this 
conclusion occurs at the cool end of the main sequence: as shown in Figure 5a, 

the is some evidence that the mean x-ray emission levels of dwarf M stars 
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drops as one moves beyond approximately dM5 (Golub 1983a), a phenomenon which 
has its counterpart in H alpha observations of such stars (Giampapa 1983). 
This suRgestion has been recently investigated by Bookbinder et al . (1984), 
who have shown for a volume-limited sample of dM stars that such a decrease in 
emission levels does indeed occur. 


5 

Rucinski (1984) has also recently considered this problem, and has come to a 

seemingly opposite conclusion (i.e., that the luminosity ratio L /l does 

X bol 

not significantly decrease within the spectrai type range M0-M6). In fact, 
there is no contradiction. Rucinski based his analysis on the sample of stars 
considered by Johnson (1983b); indeed, his Figure 1 closely resembles Figure 
5a shown here. As just discussed, the entire effect Bookbinder et al . (1984) 

find is due to the severe depletion of both detections and upper limits for 
spectral types later than roughly M5. Thus, Rucinski ‘s analysis is perfectly 
valid for the d etected stars; and his disregard of upper limits and of the 
observational decrease in the discovery rate (in either detection or upper 
bound on the x-ray emission level) of very low-mass stars makes his analysis 
not relevant to the larger question of whether late-type stars as a class 
(i.e., not just restricted to the known emitters) experience a decrease in 
activity level al roughly spectral type M5. 


The argument runs as follows. Figure 5a does not demonstrate the 
decrease in x-ray emission levels directly because of possible (and likely) 
selection effects; however, if we assume that the mass function for late-type 
stars remains roughly constant beyond M5 (Bahcall and Soneira 1980), then this 
figure does show that the x-ray observations are "misring" stars later than 
M5 . That is, since this Figure shows all the data based on beating optical 
proper motion survey data against the x-ray observations (and thus includes 
detections as well as upper bounds), and since the sum total of detections and 
upper bounds for stars later than M5 is below expectations based on a constant 
mass function, it must be that this comparision of the x-ray and optical data 
has missed a substantial number of stars. Because the residual number of 
unidentified x-ray detections in the Einst ein fields cannot make up this 
deficit, one is forced x.o conclude that either the missing stars have emission 
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levels below the thresholds of the typical Einstein fields, or the stellar 
mass function drops steeply at approximately M5V. The latter possibility 
cannot be ruled out, but in any case, the upper bound on the luminosity of the 
"missing" stars leads to the inevitable conclusion that the stars later than 
approximately M5 (if they exist in any numbers) must have 
substantially-reduced x-ray emission levels. 

{ d ) Stellar x-ray emiss i on a nd eff ective gravity. 

A particularly striking aspect of Figure 3 is that stars of all 

luminosity classes fall on the same regression line of x-ray luminosity versus 

rotation rate, indicating that stars which rotate and have outer convection 
zones emit x-rays at levels independent of effective gravity (the only major 
exceptions being the cool giants and supergiants, to be discussed below). In 
contrast, the x-ray data showed from the first observations that for early 
spectral type stars, effective gravity does play a role: that is, because 

effective gravity correlates with bolometric luminosity (and as the bolometric 
and x-ray luminosities are observational ly correlated), the effective gravity 
is perforce correlated with the x-ray luminosity. However, since theory does 
expect a relation between bolometric luminosity and x-ray emission levels (at 
least for those models which invoke instabilities in radiation-driven winds: 

see recent review articles in Underhill 1984), the correlation between x-ray 

luminosity and effective gravity for these stars is best regarded as 
secondary, and (in and of itself) of no physical interest. Unfortunately, a 
more detailed theoretical account of the x-ray emission from such stars has as 
yet proved elusive, and remains a significant task for the future. 

The theoretical situation is hardly better as one moves toward cooler 
spectral type stars. The early Eins te in Observatory data already showed that 
late spectral type giants and supergiants are at best very weak x-ray sources 
(Ayres and Linsky 1980; Vaiana 1981a, b: Vaiana et a l . 1981: Haisch and Simon 
1982). This dearth of x-ray emission from evolved cool stars remains as an 
outstanding puzzle. One popular proposal has been to invoke an analogy with 
solar coronal holes: in such regions, the coronal gas streams out to form the 

solar wind, and the level of x-ray emission is severely depressed (cf.. Maxson 
and Vaiana 1977: and the reviews in Zirker 1977). Since the evolved stars in 

question are known to be vigorous sources of low terminal speed, high mass 
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flux winds (see reviews uy Dupree 1982, Hartmann 1983, MacGregor 1983, Brown 
1984, and references therein), and as the boundary line in the H-R diagram 
defining the onset of high mass loss rates very roughly corresponds to the 
boundary line separating stars which have and have not been detected in x-rays 
(as well as in other indicators of high-temperature matter above the stellar 
surface; Linsky and Haisch 1979; Ayres and Linsky 1980; Simon, Linsky, and 
Stencel 1982), this coronal hole analogy has seemed quite apt (see also 
Cassinelli and McGregor 1984). 

The dangers in drawing this analogy too closely are severe however. 
First, it remains true that there is no commonly-accepted theoretical account 
for solar coronal holes, e.g., a first-principles theory which links a 
presumptive source of energy and momentum (such as magnetohydrodynamic waves) 
to the observed solar wind mass flux and wind speed at 1 AU, and to the 
observed temperature and density in the inner corona of the source coronal 
hole (see especially Leer, Holzer, and Fla 1982 for a general discussion of 
this problem, and Holzer. Fla. and Leer 1983 for details on the inadequacy of 
Alfven wave-driven wind models in the stellar case; see also MacGregor 1983). 
The coronal hole analogy hence is perforce purely phenomenological. Second, 
the dearth of x-ray emission from these evolved stars is far more profound 
than in the solar case: for example, the upper bounds on the stellar surface 
x-ray flux for several late spectral type supergiants lie over three orders of 
magnitude below that of solar coronal holes (Vaiana et al . 1981). This 
suggests a profound difference in the structures of the wind atmospheres of 
the Sun and these evolved stars — the steep temperature rise into the corona 
characteristic of the Sun seems to be suppressed in these far more slowly 
accelerating winds. Thus, the challenge of understanding the physics of these 
atmospheres lies In reconciling the large mass loss and low terminal speeds 
(which imply that the bulk of the flow heating and acceleration occurs 
relatively close to the star, and below the wind sonic point) with the absence 
of a temperature rise (which in turn implies either small heating rates below 
the sonic point, or very effective cooling) — a completely different agenda 
for theory than that faced by the solar physicist. Thus, to put It in its 
best light, the solar coronal hole analogy for explaining the lack of emission 
from late spectral type evolved stars remains an Interesting idea. 

( i i ) T he Pltys ica 1 C-liara ct e r i^t i cs o f S tell ar Coro na-'' . 
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Because of the essential lack of spatial resolution, it should not be 
surprising that the extent to which the physical attributes of the plasma(s) 
which produce stellar x-rays can be deduced Is severely limited. Consider 
first the coronal temperature. Even the rather modest energy resolution of 
the Einstein Observatory Imaging Proportional Counter (IPC) allowed the 

conclusion that single-component temperature models for stellar x-ray spectra 
are generally inadequate (at least for data with high signal-to-noise ) . The 
spectroscopic IPC results of the Cinsteln /Cf A stellar surveys show that one 
can grossly characterize coronal temperatures by a bimodal temperature 

distribution, with peaks in this distribution at log T ~ 6.3 and log T " 7.0 
(Vaiana 1983; Schmitt 1984; Majer et a l . 1984). Somewhat higher resolution 

spectra obtained by the Ei nstein Observatory Solid State Spectrometer (SSS) 
shows similar bimodal distributions, but with a shift in temperature: 
typically, the IPC low-temperature component is not found by the SSS analysis, 
but instead a 'mewhat high-temperature (log T > 7,5) appears (Holt et al . 

1979; Swank et al . 1981; Swank and Johnson 1982); while tne relatively high 

resolution spectroscopic data obtained from the Ei nstei n Observatory Objective 
Grating Spectrometer (OGS) for Capella (Mewe et a l . 1982) and from the EXOSAT 

OGS (Mewe 1984) seem more consistent with the IPC data. 

These results can be reconciled in the following way (Majer et al . 1984): 

any given coronal loop structure must be characterized by an internal 

temperature distribution, and hence by a fractional distribution of matter 

(the differential emission measure: Jordan 1976; Withbroe 1975) as a function 

of temperature. Furthermore, as in the Sun, it is likely that there is a 

superposition of loops with different mean properties (i.e., Ir'^gth, mean 

temperature, mean gas pressure). The result is a convolved differential 

emission measure distribution, whicfi ft>r a spatially-unresolved s teady source 

cannot be "deconvolved" into its component constituents (i.e., there is no a 

prior i method for distinguishing the contributions of the Internal run of 

temperatures within given loops from that of distinct closed structures with 

. 6 

different mean temperatures). In any case, the ultimate consequence for a 


6 

If the coronal emission varies w'th time, it may be possible to distinguish 
between different emission components (e.g., if the temporal variations have a 
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strong spectral dependence); a nice example of this possibility Is provided by 
studies of the coronal emission measure distribution on RS CVn stars (viz., 
Golub 1983a; Walter, Gibson, and Basri 1983; White et al . 1984; and the review 
by Charles 1983, and references therein). Another possibility of 
disentangling geometry already alluded to is based on correlations between 
coronal x-ray emission and emission from cooler parts of the atmosphere (viz., 
Mewe and Zwaan 1980; Ayres. Marstad, and Linsky 1981; Pallavicini et al . 1982; 
Schrijver et al . 1983; Giampapa et.^JL.. 1985), but the lack of simultaneity of 
these observations, in concert with the everpresent variability in emission 
levels, have limited use of this technique to date. 


moderate energy resolution instrument is the detection of a multi-component 
temperature distribution; quite typically, the number of such components 
corresponds to the number of well-distinguished spectral "windows" the 
spectrometer has (the IPC had two, centered at roughly 0.2 and 1.5 keV, 
respectively). Thus, the data lead to the conclusion that the coronae of 
late-type stars are characterized by a broad range of plasma temperatures, 
from log T " 6.0 to > 7.5. 

Consider next the amount of emitting material, and its density. A glance 

at the functional form of the plasma cooling function (viz., Raymond acd Smith 

- 1/2 

1977) shows that it varies roughly as T in the temperature range of 

interest to us, and with the square of the electron density; hence the cooling 
function varies at most by a factor of 3-5 due to the temperature variation 
alone. This immediately implies that the great observed variation in stellar 
x-ray luminosity is largely due to a variation in the amount of coronal 
material present . Since large-amplitude variations in x-ray luminosity 
outside of flares are rare, it is reasonable to assume that the emitting 
matter is distributed fairly homogeneously over the stellar surface. Hence 
the emitting volume is roughly fixed, and we can reasonably conclude that the 
large range of observed x-ray luminosity levels reflects a (somewhat more 
modest) variation in mean coronal density, this is precisely the case for the 
Sun (viz., Vaiana and Rosner 1978), as the modulation of coronal intensity 
within a given coronal loop structure is dictated essentially entirely by the 
variation in coronal loop density. 
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( i i i ) The Basic Physical Questions 

Can we understand the correlation between stellar x-ray emission and 
rotation for late-type stars, and between x-ray emission and bolometric 
luminosity for early-type stars? It is probably fair to say that there is now 
a consensus answer to the former question, at least on the qualitative level, 
but that a commonly-accepted model for x-ray emission from early-type stars is 
still out of reach. In the following, we briefly outline the basic model for 
stellar activity of late-type stars, and refer the reader to the 
still-developing literature on coronae of early-type stars (see especially the 
volume edited by A. B. Underhill 1985, the recent work of Owocki and Rybicki 
1984, and references therein). 

The basic notion underlying current models of stellar activity is that in 
a rotating star having an outer convection zone, coupling of rotation and 
convection with ambient magnetic fields can lead to a regenerative magnetic 
dynamo (whose ultimate driving energy source is of course the outward-flowing 
thermal flux due to nuclear burning at the stellar core; viz., Parker 1979): 
these dynamo-generated magnetic fields Inevitably rise to the surface of the 
star because of magnetic buoyancy (viz.. Acheson 1978; Schmitt and Rosner 
1983). These surface magnetic fields are still embedded in a fluid — the 
surface of the solar convection zone — in which the hydrodynamic forces 
dominate (high-beta), and hence are continuaily subjected to "jostling" by 
turbulent surface fluid motions. This "jostling" of the emerged fields can 
lead to plasma heating (for recent views, see Kuperus, lonson, and Spicer 
1981; Priest 1983; Parker 1983a, b), with the result that a corona confined by 
the emerged magnetic fields is formed; the emission from this corona is 
presumably seen in the Einstei n and EXOSAT x-ray images. 

To what extent can one improve upon this qualitative description? A 
truely first-principles theory, which starts with the dynamical equations for 
stellar structure and magnetic flux generation and dissipation, and attempts 
to predict x-ray emission levels as a function of, say, stellar mass. age. and 
rotation rate, is at present out of reach. One obvious stumbling block is our 
still meager understanding of non-linear dynamo theory, which does not allow 
us to predict with any degree of certainty the rate of magnetic flux 
production on stars, including the Sun (cf., Gilman 1983; Schussler 1983). 
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Nevertheless , there does exist a vast body of solar observations from which 
properties of solar-type dynamos may be deduced (or at least constrained). 
That is, although the dynamo itself is not directly accessible to 
observations, we can observe samples of the magnetic fields that are produced 
in the Interior, and are brought to the surface by the joint action of 
magnetic buoyancy and convection. The question is then whether the observable 
characteristics of the emerging fields show regularities which may constrain 
theory. Such studies have been conducted (cf., Golub et al . 1980; Zwaan 1980; 
Howard and LaBonte 1983; Giampapa and Rosner 1984; Vaughan 1984), and suggest 
the following picture of a solar-type dynamo. 

There are abundant reasons for believing that the dynamo is located 
rather near to the bottom of the solar convection zone, either in the bottom 
layers of the convection zone itself (Parker 1979), or in a thin 
stably-stratified undershoot layer separating the convection zone proper from 
the radiative core (Rosner 1979; Spiegel and Weiss 1980, 1982; Golub et al . 
1981; van Ballegooijan 1982; Schmitt and Rosner 1983; Schmitt, Rosner, and 
Bohn 1984). These reasons include the short time scale for buoyant rise of 
flux throughout most of the convection zone (too short for efficient dynamo 
action to proceed; Parker 1979; Acheson 1978; 1979); the effect of topological 
pumping in pushing magnetic flux downwards in the turbulent convection zone 
(Drobyshevski and Yuferev 1974; Arter 1983); the increased efficiency of 
dynamo action within a boundary layer separating regions of stable and 
unstable entropy stratification (Bullard and Gubblns 1977); and the 
possibility of stabilization of magnetic buoyancy modes in such a 
convectlvely-stable layer (Schmitt and Rosner 1983). If this is an appropriate 
model for the dynamo, there are a number of immediate deductions which follow 
when extrapolating to other stars; 

(1) For stars which are expected to be fully-convectlve on the basis of 
standard stellar interiors models, the boundary layer necessary for dynamo 
action would not exist; naively, one would thus predict that stellar activity 
ought to decrease dramatically at the point at which stars become 
ful ly-convective on the main sequence (although in the presence ol interior 
magnetic fields, it is not clear that such a star would ever become fully 
convective), indeed, a significant reduction in stellar activity levels is 
observed on the main sequence at approximately M5, the point at which such 
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stars are expected to become fully convective (but see discussion of this 
point above). 

(2) For stars with very shallow convection zones, the major effect Is 
likely to be a significant reduction in the typical size scale of surface 
activity regions (Giampapa and Rosner 1985); this result is a consequence of 
the shift to larger wavenumbers of the most unstable magnetic buoyancy modes 
(Schmitt and Rosner 1983). Again, observations (in the visible, In the UV 
continuum, and in x-rays) show that roughly blueward of the spectral type at 
which the steep reduction in convection zone depth occurs on the main 
sequence, both activity itself ( Bohm-Vitense and Dettmann 1980; Schmitt, 
Golub, Hamden, and Rosner 1985), and the amplitude modulation of activity (a 
measure of the spatial scales of surface inhomogeneities; see Giampapa and 
Rosner 1984 and references therein) show a steep decline. 

(3) For stars of fixed spectral type, an Increased gradient in 
differential rotation at the base of the convection zone ousht to have the 
consequence of significant enhancement in activity levels. As just discussed 
above, precisely such an effect appears to be present when comparing the 
activity levels of G and K stars in the Pleiades open cluster (Micela et al . 
1984, 1985). 

(4) One of the basic expectations of classical dynamo theory is that the 
efficiency of magnetic flux production is a sensitive function of the ratio of 
the stellar rotation period to the convection dverturning time scale (Durney 
and Latour 1978); this ratio (the Rossby number) enters physically because it 
is the Coriolis force that ultimately is responsible for the production of 
meridional magnetic flux from the emerging toroidal magnetic flux. Indeed, 
Noyes et al . (1984) have shown that there exists an excellent correlation 
between Ca II H and K emission strength and the Rossby number for solar-type 
stars (see also Mangenay and Praderie 1984 and Marllll and Catalano 1984). A 
more exacting test (because of the wide range of depths of the outer 
convection zone) is to consider the rotation-activity correlation for main 
sequence stars in the spectral type range A0-F7 . Thus, Schmitt, Golub, 
Narnden, and Rosner (1985) have shown, first, that the correlation between 
rotation rate and x-ray emission levels is virtually absent when either this 
full spectral type range, or subranges, are considered; but. second, that 
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wlthin the restricted subrange of F stars, there Is some evidence (as yet 
weak) for a correlation between Kossby number and activity levels. This 
seeming contradiction (which needs to be further explored by the coming 
generation of more sensitive x-ray telescopes) is precisely what is to be 
expected from a population of stars with a wide range in both rotation rates 
and convection zone depths, but for which a "rotation-activity connection" 
does hold (because the large variation in convection zone depths would mask 
any correlation between rotation rate and activity level). 


VI. STELLAR X-RAY EMISS ION I N THE MO RE GENERAL ASTROPHYSICAL CONTEXT 

In addition to providing information on the physics of the outer layers 
of stars, stellar x-ray emission — and stellar activity in general — also 
plays significant roles in the larger galactic context, some of the more 
important of which we now very briefly discuss. 

(a ) X-ray Emission from Stars and the Diffuse Soft X-ray Background 

The most straightforward consequence of stellar x-ray emission is its 

passive role in contributing to the galactic component of the diffuse soft 

x-ray background (Gorenstein and Tucker 1976; Rosner et al . 1981). This 

effect is demonstrated in Figure 6 (taken from Bockbinder et al . 1984), which 

shows the contribution of very late spectral type dwarf stars to the number 

density of soft x-ray sources as a function of detection threshold, and the 

derived contribution to the soft x-ray background (as deduced by using a 

spatial distribution for such stars derived by Bahcall and Soneira 1980). 

Comparison of Figure 6 with the results of Rosner et al . (1981), who used a 

simple exponential model for the stellar distribution, shows that the main 

conclusions remain intact; stars seem to provide approximately 20% of the 

diffuse galactic component of the soft x-ray background in the " 1/4 keV 

energy range. These results also show why it is essential that the next 

generation of very sensitive x-ray telescopes must have high spatial 

““16 *“ 2 *" 1 

resolution: at a sensitivity tlireshoJd level of 10 ergs cm s within the 
the Einstein Observatory energy passband, we expect to see about 100 dM stars 


per square degree I 
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lb) X-ray e mi ssi on From .Very Stars and the Local Interstella r 

M edium 

In addition to the above "passive" role played by stellar x-ray emission, 
it may also play an "active" role in significantly effecting the physics of 
the Interstellar medium in the immediate vicinity of stars. This possibility 
has been explored most prominently in the case of very young x-ray-active 
objects, such as are found in star formation regions: the possible effects 

which have been only tentatively explored to date include the regulation of 
star formation by the stellar x-ray emission (as a result of the large flux of 
ionizing photons emanating from these objects embedded in molecular cloud 
regions: see Lepp and McCray 1983; Krollk and Kallman 1983; Silk and Norman 

1983); the resulting increased coupling of the accreting gas to ambient 
magnetic fields (Silk and Norman 1983; Shu and Terebey 1984) implies that de- 
coupling of the accreting flow from the retarding cloud magnetic fields is no 
longer determined by the thermal history of the accreting gas viewed in 
isolation (viz., Shu and Terebey 1984, and references therein). In addition, 
these observations raise the question as to how important stellar "activity" 
was in the formative stage of the solar system and during the formative stages 
of planetary atmospheres (cf., Sekiya, Hayashi, and Nakazawa 1981: Canuto e^ 

1982: Feigelson 1982: Zahnle and Walker 1982). 

( c ) M ass and Far t i cl e Input to the Inte rst e llar Med ium 

Activity at the surface of stars is known to lead to at times vigorous 
ejection of matter from the gravitational confines of stars, either in the 
form of massive wind outflows (see reviews by Cassinelli 1979, Cassinelli and 
McGregor 1984) or transients (which themselves can lead to bulk mass ejection 
and relativistic particle production: viz., Coleman and Worden 1976) As a 

consequence of these particle and mass ejections, stars provides one of the 
main sources of mass input to the inteistel lar medium (viz. Abbott 1982 and 
Van Buren 1984), and may provide at least some of the input source particle 
spectrum for cosmic rays (either directly via the injection of particles 
accelerated during stellar flares (viz., Muilun 1979: see also Gorenstein 

1981), or indirectly via particle acceleration at shocks lying either within 
the winds themselves, or at the termination shock of the winds (viz., Cesarsky 
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and Montmerle 1983,. The latter possibilities are strongly supported by 
evidence that the elemental composition of galactic cosmic rays heavier than 
He is very similar to that of solar energetic particles (Meyer 1981, 1985); a 
recent summary of the current situation Is provided by Montmerle (1984). 


VII. CO N’CLUSIO NS AND FUTUR E P HOSPRCTS 

Our understanding of spatially unresolved astrophysical systems suffers 
from the twin difficulties that the source geometry is not obvious, and that 
the relevant physical description is not necessarily unique; these two 
problems lie at the heart of the difference between astrophysics and the 
terrestrially-based physical sciences. From this perspectice, stellar x-ray 
astronomy (at least as regards late spectral type stars) has the great 
advantage that solar observations can be used both to constrain the physics 
and to suggest the appropriate observational tools: the Sun and its environs 
provide us with a directly observable laboratory for studying the 
magnetohydrodynamics and plasma physics which must enter into the physics of 
stellar x-ray emission, and allow us to define the theoretical framework 
within which we Interpret the observations. The major results which have been 
gained from such studies are new constraints on models for stellar activity, 
which demand tests of theoretically-predicted dependences of activity levels 
with stellar parameters which cannot be varied by studying the Sun alone. 
Thus, stellar x-ray observations, and the attendant modeling, have already 
provided additional constraints on both theories of coronal heating and 
magnetic flux generation (i.e., dynamo theory) which cannot be obtained 
Independently from solar work. 

In the larger astrophysical context, observations of stellar x-ray 
emission allow us to more sensibly extend the solar analogy to more "exotic" 
(and certainly astrophys leal ly Interesting) phenomena, such as the formation 
of coronae on accretion disks and the energetics of hot galactic halos. Such 
extensions are not restricted to studies of late-type stars: certainly, an 
understanding of possible Instabilities in radiation-driven winds of OR stars 
will provide the necessary underpinnings for any attempts to appiy such 
theories to, for example, instabilities in radiation-driven jets. 


-34- 


Where does the future of x-ray observations of stars lie? A strong hint 
of the Most fruitful lines of research is easily gleaned by considering the 
Major outstanding observational probleMS. each of which is intimately tied 
a major theoretical problem of stellar activity: 

{ a ) What Ij the nature of the emission from early spectral type stars ? 

That is, at present neither the true source spectrum, nor the nature of the 

low energy cutoff seen in low and moderate resolution spectra, are 

well-understood. At issue here is the physical location of the x-ray-emitting 

matter (especially relative to the stellar wind outflow), and whether the 

observed emission is best regarded as a remote analogue of solar magnetic 

surface activity — a corona -- at or very near the stellar surface (Hearn 

1972: Cssslnelll and Olson 1979; Rosner and Valana 1980; Stewart and Fabian 

1981; Cassinelli 1984; Uchlda 1984; Waldron 1984; Underhill 1984), or as a 

consequence of instabilities in the massive radiation-driven winds (Lucy and 

Solomon 1970; Hearn 1975, 1984: Nelson and Hearn 1978; McGregor, Hartmann, and 

Raymond 1979; Carloerg 1980; Lucy and White 1980; Lucy 1982; Owocki and 

Rybicki 1984). These problems will be directly addressed as multi-wavelength 

observations of the nonthcrmal activity of such stars proceed (cf., Caillault 

et al . 1985; White 1985; and references therein), and at x-ray wavelengths 

will require the kind of very high sensitivity, high resolution spectroscopy 

which currently-planned instruments such as AXAP will permit. High 

sensitivity is essential because the nearest subject stars He roughly 400 pc 

30 32 -1 

from us, and have x-ray luminosities of order '10 - 10 ergs sec ; this 

3 

implies photon-limited observations for the kind of high resolution (> 10 ) 
spectroscopy needed to address the above problems. 

(b) Do very l ate sp ectral type gian ts and supergiants emit any x- rays? 
The lack of x-ray detections of such stars suggest that either they do not 
emit x-rays (or do so only at levels substantially below solar values), or 
that extinction (by, for example, stellar wind material) is so prohibitive 
that there is no significant emergent x-ray flux. The basic question here is 
whether the solar analogy with coronal holes (regions of sharply-reduced x-ray 
emission which are the sources of the solar hig)i speed wind streams) is 
appropriate, so that the lack of observed emission accurately reflects an 
absence of any x-ray emitting matter (Ayres and Linsky 1980; Dupree 1982; 
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Linsky 1984). The future observational requirements are thus to better 
constrain the x-ray spectra of stars near the "dividing line" (Linsky and 
Haisch 1979) in the H-R diagram (in order to test for a softening of the 
spectra upon reaching the dividing line, and to search for a soft cut-off 
which may result from the presence of intervening cool wind plasma); and to 
obtain far stronger upper limits on the possible emission levels of late 
spectral type giants and supergiants. Translated into lne.crument.il terms, the 
observational needs are again for high resolution, high sensitivity 
spectroscopy and for very high sensitivity, high spatial resolution imaging. 

( c ) Do normal main se que nce stars in the spectral type range 

approximately B5-A4 emi t an^ x-r ays ? Present observations show that x-ray 

levels of such stars (if they emit at all) must lie substantially below the 

median level of solar-type stars, and that a steep rise in emission levels 

occurs in parallel with the sharp growth in depth of the outer convection zone 

as one moves down the main sequence past spectral type approximately A4 (a 

boundary line whose precise location predicted by stellar Interiors models is 

sensitive to the value of the mixing length parameter — the ratio of the 

mixing length to the local scale height). However, the upper bounds on x-ray 

emission levels for late B and early A dwarf stars do not lie very much below 

solar emission levels, and it is of great interest to place much tighter 

constraints on the possible emission levels: we want to know to what extent 

stellar activity has as a necessary prerequisite a convective stellar surface 

(which the stars in question are presumed not to have). The key observational 

-16 

requirement in this case is clearly sensitivity; at a limiting flux of 10 
-1 -2 

ergs sec cm (which will be easily accomplished by ^ A£) . future 
observations will be able to place an upper limit on the coronal x-ray 
luminosity of the nearest candidate stars of less than log L " 23 5 more 

X 

than a factor of 1,000 below the typical solar x-ray luminosity. 

( d ) What is the nat ur e of the x-ra y emission from young (pre-main 

sequence) stars ? Thus, as is the case for the early spectral type stars, at 
present we know neittier what the true source spectrum is, nor what the nature 
of the low energy cutoff seen in low resolution spectra is. Furthermore, the 
nature of the observed strong variability is highly uncertain (althjugh an 
argument based on the flare energy-frequency histogram suggests a solar flare 
analogy; Montmerle 1983; Rosner and \aiana 1978). Here the central 
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problems are to what extent the solar analogy applies; whether there is any 
quiescent component; whether the emission is best seen as (solar) flare-like; 
and to what extent the observed "activity" is the result of a (regenerative) 
dynamo process, as opposed to the decay of primordial fields. Again, high 
resolution spectroscopy appear essential to answering the question of the 
structure of the emitting regions. Furthermore, because these objects 

typically lie in regions of large obscuration at soft x-ray wavelengths, and 
as their spectra are somewhat harder than those of main sequence stars of 
comparable spectral type, there is considerable reason to study these objects 
with imaging telescopes capable of high sensitivity at energies above 2-3 keV 
(imaging is essential because the projected number density of objects in 
star-forming regions can lead to severe source confusion at angular 

resolutions as good as 1'). 

( e ) Is the dearth of x-ray-brigh t (a ctive) M dwarf stars later than 

approximately M5 a reflection of a fundamental change in dynamo behavior for 
fully convective stars, or is it simply that these "missing" stars are truely 
missing? If the latter is the case, then the mass function at the low-mass 
end of the main sequence must steeply decline past roughly M5, with major 
implications for either theories of low-mass star formation, or for models of 
stellar evolution (which commonly predict a very slow descent to the main 
sequence — and hydrogen burning -- for these stars). From the x-ray point of 
view, the key task is to probe the solar neighborhood for very low-luminosity 
stellar objects; if the "missing" stars do exist, then one would expect to 
find them in such "deep" surveys. The key observational requirements are for 
high sensitivity and high angular resolution imaging; the latter requirement 
is borne of the fact that the optical counterparts to the sought-for x-ray 
sources will be very faint stars (ra > 22), so that source identification will 

V 

be extremely difficult unless the source position error circle is sufficiently 
small (i.e., of order an arc second, or less). 

( f ) Is the transient component of stellar x-ray emission a reflection o f 

a fundamentally dis t inc t plasm a heating process? As discussed above, 

transient emission appears to dominate the x-ray output of some classes of 
stars, and in fact constitutes the entire emission from pre-main sequence 
stars observed to date. Because the frequency distributions of flares 
generally show a power law dependence on total flare energy or peak flux for 
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many classes of transient sources (viz., Rosner and Vaiana 1978), it is not at 
all obvious whether such transients contribute significantly to what is 
normally referred to as the quiascent emission component: the crucial question 
is then whether there is a iower energy cutoff on this power law. There is 
some Indication from observations of solar hard x-ray emission that such a 
cutoff (if it exists) may lie at sufficiently low flare energies that the 
total energy associated with fast electrons produced in "microflares" is a 
significant fraction of the total coronal x-ray luminosity (viz., Lin et al. 
1984). If this is the case, then the traditional distinction between 
transient and quiescent energy release and plasma heating disappears, and we 
will have to regard stellar flares as simply the extreme high energy tall of a 
continuous distribution of coronal energy release events. Whether this is in 
fact the case remains to be established, and it is clearly solar observations 
that will provide the definitive resolution to this problem. 


These few examples should make it plain that it is indeed both 
spectroscopy and very high sensitivity and spatial resolution imaging which 
are essential to the construction of activity models more complex than the 
simple models discussed to date. It is likely that the next generation of 
x-ray telescopes (such as a modified ROSAT . or NASA-sponsored high-throughput 
missions) and facilities (i.e., AXAF) will indeed carry out the necessary 
observations . 
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LI C URE CAPTIONS 


Figure 1 : Schematized indication of where stellar x-ray emission has been 

detected in the Hertzsprung-Russell diagram. 

Figure 2: Integral x-ray luminosity functions for F, G, and M stars on the 

main sequence, as derived from volume-licited samples obierved by the 
Einstein Observatory (Rosner et al . 1981; Bookbinder, Golub, Schmitt, and 
Rosner 1984; Maggio et al ■ 1984; Majer et al . 1984; Schmitt et al . 1984); 
the Sun's coronal luminosity places It near the bottom of this 
distribution of stars in this spectral type range. 

Figure 3: Scatter diagram of soft x-ray luminosity versus rotation rate, 

adapted from Pallavicini et al . (1981, 1982), and further updated using 

both newly-analyzed Einstein data as well as recent EXOSAT data (courtesy 
J. Schmitt). O nly st ars with either measured or bounded rotation rates 
are show n. Note that in spite of the diversity of stellar types in this 
sample, the correlation between x-ray luminosity and rotation period 
remains very strong. 

Fig u re 4 : Cumulative x-ray luminosity functions for G stars in the Pleiades 

and Hyades open clusters, and for dwarf G field stars; note the clear age 
effect — the cumulative x-ray luminosity function shows a systematic 
shift to lower emission levels as the stars age. 

Fig u re 5: X-ray emission at the red end of the main sequence; (a) scatter 

diagram for the variation of x-ray luminosity with spectral type in the 
interval M0-M8; (b) comparison of the derived cumulative x-ray luminosity 
functions for M stars < M5 and > M5. See text for discussion. 

Figure 6: In this figure, we plot (solid curve) the log of the predicted 

number of stars observed per square degree whose x-ray flux lies above a 
threshold value S, log N(>S), versus the log of this threshold flux, as 
derived for stars in the spectral type range M0-M5 from the pointed 
Einstein Observatory stellar data (i.e,, using the luminosity function 
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shown in Fisure 6. and the M dwarf spatial distribution in our (raiaxy 
from the model of Bahcali and Soneira 1980). Also shown are (a) the 
observed values of log N(>S)-log S for all stars derived from the 
Einste in Observatory Medium Sensitivity Survey (Maccacaro et al . 1982; 
Stocke et al . 1983), and (b) the log N(>S)-log S curves (dashed curves) 
for stars in the spectral type range M5-M8 under two distinct 
assumptions: (i) their x-ray luminosity function is identical to that of 
the earlier spectral type, and (ii) the observed sample is a fair 
representative of the universe of such M dwarf stars (so that the x-ray 
luminosity function derived directly from the data is used). Note that 
the stars in the range M0-M5 already satisfy the observed x-ray log N-log 
S behavior, so that the stars in the range M5-M8 cannot possibly be 
vigorous x-ray emitters. The results shown here are taken from 
Bookbinder et al . (1984). 
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The Time Resolution Domain of Stellar Radio Astronomy 


Jay Bookbinder 
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Harvard University 


I. Introduction 

While the origins of stellar radio astronomy have their roots in time-variable 
phenomena, high time resolution (HTR) radio observations of stellar sources is a 
very youn£ technique. This is a somewhat surprising development in light of 
the historical origin of stellar radio astronomy, which was exclusively con- 
cerned with the flare activity associated with dMe stars. One would expect that, 
having detected the flare emission, higher time resolution would have been em- 
ployed to elucidate the time structure of the flare, thereby providing constraints 
on the emission and particle acceleration mechanisms. Rather, focus shifted from 
flare observations to detecting a variety of stellar sources at ever-lower flux 
levels, so that "quiescent" emission has been detected from severed areas on the 
HR diagram. In fact, we shall see below that it may be possible that high time 
resolution studies can clarify the nature of the "quiescent" emission in some 
stellar sources. 

A number of basic questions are appropriate for discussion at this point. 
What constitutes high time resolution in the context of stellar radio observa- 
tions? For what sources are these observations likely to yield interesting 
results? And what current problems in understanding the emission can best be 
addressed by these observations? Can we use the solar analogue, and in particu- 
lar the concept of the "solar-stellar connection", as a guide for our observations? 
Is it possible for these observations to provide us with a better understanding of 
related stellar properties, i.e. such as the nature of the stellar d>iiamo? Finally, 
are the facilities that are currently available or will be a' ailable in the near 
future suitable for making these observations? 

For the purposes of our discussion, high time resolution refers to timescales 
that are short compared to either the hydrodynamic flow timescale or the ther- 
mal cooling timescale. We are thus selecting emission processes that are typi- 
cally called "non-thermal", even if the underlying electron distribution is "ther- 
mal", i.e. Maxwellian. With this criteria, we eliminate from our discussion the 
free-free brtsmmstrahlung emission from, e.g., most 0 stars, which are variable, 
for example, on the timescales of months (see however, Abbott, Bieging and 
Church well 1984 and White 1984). 

One consequence of higlier time resolution is a better estimate for the 
brightness temperature, that is determined for the source. In general, one 
assumes that the timescale of variability represents seme real propagation effect 
of the exciting phenomenon, and thus a characteristic velocity. This yields an 
upper limit on the true source size, which is simply related to the brightness 
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temperature by: 


Tg ■ 6.41 X 10'° (S f/f^) (d/R)^ 

where Is the radio flux at frequency f, d is the distance to the star In pc and 
R is the source size In. units of solar radii. Obviously, low time resolution 
observations only provide a lower limit on the source's true brightness tempera- 
ture. If R^ is to be reasonably constrained for most stellar flare observations 
(using the rise-time of the flare), then the integration times must not exceed 
about one second for M dwarfs. A second effect of long integration times is to 
reduce the peak fluxes by averaging the instantaneous fluxes over the integra- 
tion time. One example of this effect is found in the observations of Gary, Lin- 
sky and Dulk (1982) of a flare on L726-A (one of the components of UV Ceti). 
Their five minutes integrations gave a peak flux of under 12 mJy, while 10 
second integrations had peak fluxes over 17 mJy and showed evidence that 
finer, but unresolved, time structure existed. The Importance of a strong con- 
straint on Tg cannot be over-em^asized in terms of limiting the possible emis- 
sion mechainsms: if Tg > 3 x 10° , the emission is prob^,bly coherent; if Tg > 
10^^, then there are no possible incoherent emission mechanisms. By coin- 
cidence, these short timescales are also near the observational limits imposed by 
the current generation of radio telescopes. In general, observations on longer 
timescales have few instrumental restrictions since long integration times do 
not hamper variability studies on those timescales. 

Given the above restrictions on the emission mechanisms and timescales, 
what stellar sources lend themselves to HTR observations? Besides the sun, Gib- 
son (1984) has compiled a list of 77 stellar sources that show evidence of 
"nou-thermal" emission. Among these, nine are listed as sources of coherent 
radiation, and another nine have been observed with both coherent and 
incoherent emission. Each of these sources is therefore a good candidate for HTR 
observations. These sources cover a wide range of the HR diagram, with some 
significant gaps, but are mostly M giants, flare stars (M dwarfs), RS CVn sys- 
tems and Wolf-Rayet and 0 stars and Algol-type binaries. We will restrict our 
discussion to the late- type stars and to the RS CVn systems, even though many 
of the points made in the foliowring dlscu.‘>sion are valid for the more exotic 
sources such as SS433, Cyg X-3 and Cir X- 1. We wrill discuss as examples some 
of these sources and identify what information HTR observations can provide. 

II. Solar-type Stars 

Flare stars have been the object of radio variability studies for over twenty 
years, since Lovell (1963) first observed stellar flares. Since then, observations 
have been made from below 20 Mhz to about 15 Ghz, with a variety of sensi- 
tivities and time resolutions. For these stars in particular, the \ise of the solar 
analogy may be quite apt, given that the overall structure of the stars is quite 
similar to the sun. One caveat, based upon solar analogy, is that in interpreting 
stellar observations we must expect short-lived radio phenomena may have 
vastly different properties, and may be sampling extremely different plasma 
conditions, depending on the wavelength of the observation. While the earliest 
detections of stellar flares occured at long wavelengths(430 MHz), most recent 
work has been done in the 1.4 Ghz to 5 GHz bands, especially since the advent 
of the VLA. 
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It is well accepted that solar and stellar activity at the chromospheric and 
coronal levels is a product of the magnetic field emerging from the stellar inte- 
rior, and that this activity (at least in the solar case) is closely correlated urith 
the active regions • i.e. spots and plages. Indirectly then, the observation and 
characterization of stellar flares may provide, via the constraints placed upon 
the magnetic fields involved with the flare, insights on the nature of stellar 
activity. One of the most pressing questions that can be asked (and answered 
only by stellar observations) is whether the near-equipartltlon of gas and mag- 
netic pressure that Is maintained by the sun is Inherent to all solar-type stars 
(I.e., having surface convection zones). If so, then how does the magnetic field 
depend in detail on the two parameters most closely associated with stellar 
dynamos - the age and rotation rate of the star? To answer these questions will 
require a large amount of observing time dedicated to the task of establishing 
the statistical distribution of stellar flare energy (and perhaps other characteris- 
tics as well) a function of event frequency from well chosen samples. Proper 
analysis of these distributions, to determine whether they can be represented by 
the solar distribution (and if not, determining which stellar parameters are 
responsible for the difference) will require the use of sophisticated statistical 
techniques. Thus, efforts should be made to extend detections of flares to all 
spectral lypes that are suspected of possessing solar-type closed magnetic struc- 
tures (loops), since flares are typically associated with these features. 

It is Interesting to compare the highest time resolution observations yet 
reported for stellar radio flares with a solar example. Lang et aJ. (1983) 
observed an extremely strong flare (peak fluxes above 120 mJy at 1.4 Ghz) on 
AD Leo that contained very fine structures that were unresolved on timescales 
of 0.2 seconds (see Figure la). Similar solar observations (Slottje 1978,1980) 
show complex structure down to the millisecond level (Figure lb) that is quite 
similar in appearance to the stellar flare. The importance of the I..ang et al. 
observation lies in the high degrees of polarization and deduced brightness tem- 
peratures in excess of 10^*^ K. the first indication that coherent processes are 
responsible for (some of) the emission at this frequency. Care must be taken in 
making comparisons to solar microwave bursts, however. While the typical 
solar burst consists of an impulsive phase followed by a gradual component, the 
gradual component that followed the extremely spiky emission on AD Leo could 
not have been the stellar analog of the solar gradual phase (cf. Holman, Book- 
binder and Golub 1984). 

Longer wavelength observations also have showu the existence of fine time 
structure in both the solar and stellar casf.s. Solar bursts with fine structure on 
the order of 1 - 100 msec are observed at a wide variety of frequencies: below 
1 Ghz (Slottje 1972), 1 to 1.4 GHz (Droge 1977), 2.6 GHz (Slottje 1978), 10 GHz 
(Hurford et al. 1979) and 22 and 44 GHz (Kaufman et aJ. 1980). Similarly, stel- 
lar bursts have been observed (usually with durations of less than a minute and 
variability on the order of a few seconds) at 196, 318 and 430 MHz (Spangler 
and Moffett 1976, Spangler, Shawhan and Rankin 1974) and 1.4 GHz (Lang et 
el. 1983). Since we are Interested In the solar analogy, it may be worth noting 
from a solar example that variability may be best found by studying the polar- 
ized emission - in one example of a Type III b’.TTst, essentially no detailed struc- 
ture was evident in the total intensity, but was present in the polarization (Pick 
et el. 1960, also Slottje 1974). 

Furthermore, the timescale for the variability, if it does represent a real 
propagation effect (such as Alfven waves or slow-mode MIID waves), is 
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interesting for another reason as well, since It then measures the slr^ of the 
coronal structure ^vln^ rise to the emission. The source dimensions can be 
estimated by assuming a spatial ^owth rate between the assumed propagation 
speed (e.g. the Alfven speed) and c. These would be the first measurements 
(albeit indirect) of the sizes of coronal loops on stars other than the sun. 
Through the use of HTR obervatlons, even without detailed knowldege of the 
source spectra and polarization, one can learn something about the microphysics 
of the flare plasma - l.e., the electron density, the magnetic field, the particle 
energy, etc (cf. Holman, Bookbinder and Giolub 1984). The use of these diagnos- 
tics on a variety of active late-type stars, from dM through F (all of which are 
suspected to possess solar-llke uiosed magnetic structures) can provide clues on 
the effects of the stellar parameters (such as efectlve temperature and surface 
gravity) on the stellar loops. 
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Fig. la. (above) A recording of the left circular polari- 
zation of a large flare observed on AD Leo at 1.4 GHz 
with 200 msec integration times. Note the extensive 
(and unresolved) structure (from Lang et al. 1983). The 
R.C.P showed no such variability, much like Its coun- 
terpart In Fig. lb. 

Fig. Ib.(left) A solar observation at 2.65 GHz (after 
Slottje 1980). The two figures are displayed with the 
same scales for their abclssas. The similarities are strik- 
ing. 


Another Important consequence of HTR observations are the restrictions on 
the particle acceleration mechanisms: the rise and decay times of the emission 
and (In the case of spike) the Intersplke period provide Information. Most 
models prior to Slottje (1976) did not require acceleration mechanisms to be 
affective on timescales less than about a second. A number of theoretical 
mechanisms are available that can provide hi^ brightness temperatures and 
.spiky emission, such as ^yrosynchrotron maserlnA. stimulated plasma emission 
from the coherent Interaction of electrostatic upper hybrid waves and the 
Interaction of electron plasma waves (Holman 1983 and references therein). If 
we consider an electron-cyclotron maser, for example, the rise time of the emis- 
sion reflects either the Increasing number of electrons In the trap (l.e. the dura- 
tion of the electron pulse), the saturation time of the maser, or else simply vari- 
ations In the electric field Itself. The decay time reflects either the decreasing 
heating rate or the e-folding time for the electrons to precipitate, whichever Is 
shorter. The Inter spike period then Is a measure of the timescale required to 
accelerate the electrons, the loop transit time or, more likely, represents the 
time during which the rather stringent conditions required for maser operation 
are not satisfied. Moreover, accurate knowledge of the brightness temperature 
can also specify the temporal form of the maser. For a maser with emission at 
the second harmonic to operate In a steady state condition requires that Tg < 

K, though an Intermittent maser can produce spikes with T^ up to 10^ K 
(Melrose and Dulk 1962). 

One of the principal difficulties In studying the characteristics of stellar 
flares Is the relatively low rate of strong flares from these stars - and of course, 
as one moves to smaller flares, one runs Into the additional problem of dlstln- 
guisbing a flare from the noise in the signal. Fisher and Gibson (1981) provide 
rough estimates of flare frequencies (at 1.4 and 5 GHz) on the stars in the UV 
Cetl system: roughly one per 3.5 hours of observing. I have estimated the rate 
for AD Leo as one flare (greater than 20 mJy) per 8 hours at 1.4 GHz. LoveU 
(1963) has reported a flare rate of one flare per 35 hours at 240 Mhz, though 
this result is for very strong flares (about 1 Jy). From Spangler and Moffett 
(1976) and Spangler, Shawhan and Rankin (1974) we can estimate stellar flare 
frequencies at longer wavelengths for some of the more active M dwarfs For 
Wolf 424 and YZ CMl, we obtain mean times between flares (greater than 0.5 
Jy at 196 and 318 MHz and 1.5 Jy at 430 MHz) of 12,2.4,1 hours and 
2. 8, 1.8, 3 hours at 196,318 and 430 MHz respectively. 


111. HS CVn Systems 

Having looked at the information HTR can provide for us with respect to the 
solar-type stars, what can be said about the future of HTR observations for RS 
CVn systems? First, we note tnat HTR as applied to these objects involves 
timescales significantly longer than for the dM stars - and hence much easier to 
achieve In practice. Thotigh the flare events generally occur on the timescales 
of hours to days, we note that observations of shorter-lived phemomeuon are in 
the literature. Brown and Crane (1978) observed variations in the circularly 
polarized 2695 Mhz flux from V711 Tau on the order of a few minutes. Fix et 
al. (1980) have also observed variations of factors of two In the 1665 MHz 
flux on the order of several minutes. These observations provide us with an 
example of the different temporal behavior that can be found at different fre- 
quencies; Feldman (1983) reports that over 1000 hours of observing RS CVn 


systems dt 3 cm has sho(vu none of the rapidly changing (highly) circularly 
polarized fluxes that have been reported at longer wavelengths (?? cm and ?? 
cm, respectively) by Feldman (1978) and HJellming and Gibson (1980). The Y 
implication is that Incoherent gyrosynchrotron emission seems adequate to 
explain the 3 cm emission, but that there is evidence for more exotic (coherent 
gyrosynchrotron or plasma emission) processes at the longer wavelengths. 

It is also possible that HTR can be used to establish the nature of the quies- 
cent emission from these systems. Mullan (1984) has indicated that, unlike for 
the flare stars, the flare-frequency vs. flare-energy curve of RS CVn systems is 
consistent with the hypothesis that their "quiescent" emission is really the 
superposition of a large number of small flares. More high time resolution 
observations could test tb<s hypothesis more rigorously by searching for evi- 
dence of large numbers of small flares. HTR observations of these systems on 
timescales much less than a second will likely be useful only to study the ini- 
tial phase of the larger flares, when the primary energy release is occurring, 
unless the flares really are a manifestation of the overall coronal heating pro- 
cess. Along the same lines, but at the other end of the energy scale, are the so- 
called "superflares", such as reported by Simon et al. (1980) Feldman (1963) and 
HJellming and Gibson (1980), with flux densities on the order of BOO millijan- 
skys (at 10.5 Ghz), which are prime candidates for HTR observations. Fine 
structure, on the order of tens of milliseconds, should be resolvable il it exists. 
Although the smoothness of the rise and fall of the so-called superflares argues 
strongly against an intermittently operating mechanism, higher time resolution 
of these flares might provide interesting counterpoint to the current theory. 

There also exists one other controversial question about the RS CVn systems 
that can be addressed by these same observations - whether the cool spot that 
occurs represents a single active region or an association o^ several smaller, but 
energetically "independent" regions. Rosner and Vaiana (1978) have shown that 
it is possible to fit power laws to the flare-energy vs. flare-frequency curves 
for flares from many active stars under the assumption that the flaring is a sto- 
chastic relaxation phenomenon. Their model allows one to make inferences on 
the mode of energy storage and release. Such a fit to the radio flares from an 
RS CVn would indicate that the region is a single unit form the point of view 
of its energetics. However, if it is not possible to fit a single power law to the 
curve, then the spot is an agglomeration of several energetically Independent 
active regions (i.e. the total energy involved in a flare should be related to the 
Interflare duration). The use of this 2 uialysis is also pertinent to the question of 
the quiescent emission as discussed above. 


V. iDstrumenlatiou 

As we have pointed out above, the high time resolution that is needed in stellar 
obseivations must compete against the need to have sufficient signal to noise 
(i.e. long integration times). In general, stellar observations should be made at 
the highest possible time resolution and then Intergrated up to provide suffi- 
cient signal-to-noise. Some single dish instruments (i.e. the 305m dish at Are- 
cibo) can succesfully reduce their integration times to the region of interest (on 
the order of a tenth of a second or less), provided that signal has a reasonably 
high peak flux (cf. liaug et al. 1983, Bolce et al. 1981). As usual with single 
dishes, confusion and terrestrial interference might be considered problems. 
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However, in searchiug for variability at these timescales the problem of confus- 
ing sources is negligible, and interference monitors have proven quite effective 
in discriminating between local signals and those of stellar origin (Spangler and 
Moffett 1976). Recently, proposals have been made (Stinebring 1984) that will 
permit extremely high time resolution observations using the VLA as a phased 
array, reaching down to about 20 microseconds from the current 3.3 second 
minimum integration time. Use of the VLA in this mode should be encouraged 
for all stellar observations, since it is run parallel *o the aperture synthesis 
mode; i.e. maps can be made simultaneously with high time resolution flux 
measurements. It is worth stressing that the instrumental requirements for 
HTR radio observations do not demand major technological advances and can be 
easily met by existing technology (though perhaps at substantial costs). 


VI. Summary 

In conclusion then, we find that one of the main limitations on HTR observa- 
tions is the refusal of stars to cooperate with the observer - i.e. low flare rates. 
Instrumental problems, mostly obtaining the necessary sensitivity (i.e. low 
noise) on short integration times is also a major problem, though not technically 
insurmountable. Few Instruments are likely to be able to excel Axecibo or the 
VLA for this mode of observing. High time resolution observations are neces- 
sary to determine the nature of both the acceleration and emission mechanisms 
responsible for the short-lived radio phenomena that have already been 
observed. Further observations will undoubtedly provide us with novel events 
not familiar to us from the solar context. Besides the usual flare stcu-s and RS 
CVn systems, it may be of interest to perform HTR observations on objects such 
as W UMa stars, PMS objects, OB stars, AM Her systems etc., all of which show 
evidence for non-theimal emission and/or short timescale variability. 
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ON MAGNETIC FIELD STOCHA5TICITY AND NON-THERMAL UNE BROADENING 

IN SOLAR FLARES 
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ABSTRACT 

We discuss observations of non-thermal line broadening seen in solar flares by the Solar 
Maximum Mission Satellite In light of recent results on the generation of magnetic field 
stochasticity. We show that a consistent model for the data can be constructed by as- 
suming that the observations signal the destruction of an ambient magnetohydrodynamlc 
equilibrium. 
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I. INTRODUCTION 

The spectroscopic capabilities of the 5oiar Maximum Mission (SMM), P78-1, and 
Hinotori satellites have led to a number of qualitatively new observations of solar 
flares, perhaps the most striking of which are the data on line broadening observed dur- 
ing the initial stages of flares (viz., Culhane et al. 1981, Tanaka et al. 1982, Doschek et 
al. 1985). In this paper we shall focus on one particular type of line broadening obser- 
vation which we believe may hold a key to the processes which lead to the disruptive 
behavior characteristic of solar flares. 

The onset of solar flares is characterized by large nn»-thermal broadening of soft 
jc-ray lines emitted by highly ionized heavy ions, such as Ca XIX and Fe XXV. This 
effect is usually attributed to the presence of non-thermal plasma motions, and may be 
observed one or two minutes before the impulsive increase in hard x-ray flux (which 
traditionally marks the onset of flaring), and before the appearance of high-speed 
upflows of chromospheric material heated to corona temperatures (chro^-iospheric "eva- 
poration"; Antonucci et al. 1982, 1984a). The excess line widths persist as long as 
there is observational evidence for energy release in the flare site. In the decay phase 
of flares, non-thermal velocities either are not observed, or are present at very low lev- 
els. The turbulence level in the plasma appears to be independent of the position of the 
flare on the solar disk. 

A number of Interpretations of these data are possible, the most commonly-accepted 
of which regard the overall line profile as a superposition of r'^arious Doppler-shifted 
components; this superposition is believed to be a result of an integration along the 
line-of-sight of various dLstinct, unresolved loop structures (or parts of loop structures). 
In this paper, we wish to pursue an alternative picture, in which the observed line 
broadening is due to a superposition of Doppler-shifted line profiles arising from distinct 
plasma flows originating within a single loop structure. We shall show that this 
hypothesis not only gives good account of the data, but also that it is a natural theoreti- 
cal consequence of the lack of topological stability (viz., Tsinganos, Distler, and Rosner 
1984) of coronal magnetic structures. 

The structure of our paper is thus as follows. We first present a brief summary of 
the data (^ II), describe the model and apply it to the data (^ III), and then summarize 
our conclusions ($ IV). 


II. THK DATA 

Higli resolution spectroscopic observations of the thermal plasma formed during 
solar flares have been performed during the last solar maximum with the Solflex 


expeiiineiit on the P78-1 satellite (Doschek et al. 1979). the Soft X-ray Polychromator 
on board the SMM satellite (Acton et al. 1980). and the SOX spectrometers on the Hino~ 
tori satellite (Tanaka et al. 1982). The data from these observations have shown that 
Jaat before and during the impulsive phase (when energy is presumably first 
released), the width of soft X-ray spectral lines is larger than expected from the ther- 
mal broadening, as predicted from the electron temperature derived for the emitting 
plasma. This broadening effect is seen throughout the impulsive phase; however, the 
line broadening does decrease as the impulsive phase wanes, so that lines are thermal (or 
almost thermal) after the peak in soft X-ray flux (Doschek et al. 1980, Feldman et al. 
1980, Antonucci et al. 1982, 1984a, Tanaka et al. 1982). These results are reviewed by 
Doschek et al. (1985). 

The spectral regions preferentially observed by these experiments are those of the 
Ca XIX and Fe XXV helium-like resonance lines, formed at temperatures typical of flare 
plasmas (which are in the range 1-3x10 K). The satellite lines formed in the same 

spectral regions provide direct information on the electron temperature of the source 
(Gabriel 1972). Since the electron temperature is well-determined, the non-thermal 
excess in line width can be confidently derived. This excess cannot in general be attri- 
buted to a difference in ion eud electron temperatures, became at the densities estimated 
for flare plasmas, the reqtnred temperature difference cannot be maintained for time- 
intervals as long as the flare rise times. The estimated turbulent velocity amplitude 
derived from the presumed Doppler temperature is between 100 to 200 km sec'^. Non- 
thermal broadenlngs are also observed in cooler coronal lines and in transition lines 
(Doschek et al. 1985) during the impulsive phase. 

Analysis of a large number of flares observed with the Bent Crystal Spectrometer of 
the SMM Soft X-ray Polychromator (Antonucci et al. 1984a) revealed the following key 
features: 

1. The noiz-thermal excess in line width is larger at flare onset, and 
decreases monotonically during the impulsive phase. Thus, the maximum 
broadening is found to precede or to coincide with the peak in hard X-ray flux; in 
the decay phase, spectra) lines are thermal, and late in the decay, the line width 
may increase again, with a velocity parameter of about 60 km sec"^. 

2. The degree of non-thermal excess appears to be uncorrelated with the 
position of flares on the solar disk. Thus, an average non-thermal velocity of 
100±40 km sec’^ is found for disk flares, and 120^30 km sec”^ for flares at longi- 
tudes beyond 60°. Furthermore, limb flares (such as the flare of June 29, 1980) 
show non-thermal excesses comparable to that of flares observed on the disk. 
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3. When significant soft X-ray emission is detected before the hard X-ray 
burst, the non-thermal line broadening is observed to increase one or two minutes 
before the hard X-ray burst. The increase in line width is in general accompanied 
by moderate intensity non-lmpulslve hard X-ray flux. This increase in line width 
for disk flares occurs before blue-shifted components are observed in the spectra 
(note that the blue shifts are observed to be simultaneously with the impulsive 
increase in hard X-ray emission at burst onset). Those blue shifts are usually inter- 
preted as the signature of high-speed upflows of chromospheric material as it Is 
heated to coional temperatures during the "evaporation" process (Antonucci et al. 
1982). 

The particularlly well -observed flare of May 21, 1980 tjrpifies these characteristics. 
In this event, which was not preceded by emy detectable soft X-ray emission at tempera- 

7 

tures of 10 K from the active region, it is possible to clearly distinguish the initial 
phase, characterized by non-thermal velocities of about 220 km sec~^ and moderate 
intensity hard X-ray flux, from the impulsive increase in hard X-ray emission and the 
onset of plasma upflows, as shown in figure 1. At this stage, soft X-rays are emitted 
predominantly from one of the two footpoints of the flaring loop system (Antonucci et 
al. 1984’ , c). Flare footpoints can be cl arly identified in a second stag. , in coincidence 
with the impulsive increase in hard X-ray flux, whan localized hard X-ray sources 
appear in correspondence with the ila kernels. This emission is thought to be due to 
energy deposition at the loop base by non-thermal electrons (Hoyng et al, 1981). 

The first spectrum obtained during the May 2 1 event in the Ca XIX spectral region, 
averaged over a two-minute period, is shown in Figure 2a, The initial time of the 
observation period is 20:53:02 UT. Spectral features are not resolved due to the large 
non-thermal broadeuings, and the profile of the resonance line is consistent with a 

U 

Doppler temperature of 1.3x10 K, while the inferred electron temperature is only 
9x10 K, The difference in Doppler and electron temperatures implies a non-thermal 
velocity of 220 ± 30 km sec'^. This Initial Calcium spectrum is to be compared with 
the spectrum obtained at the peak of the s«}i't X-ray flare (Figure 2b). The continuous 
curves superimposed on the data are 'ynthotic spectra computed as described in 
Antonucci et al. (1882). 

The soft X-ray image in the energy band from 3.6 to 8 keV of the flare region in 
the initial state at 20:53:02 UT (as obtained with the SMM Hard X-ray Imaging Spec- 
trometer; van Reek et al. 1980) i? shown in Figure 3. The spectrometer has a spatial 
resolution of 8". Points A and B in Figure 3 indicate the position of the footpoints of 
the flaring loop system inferred from the hard X-ray images in the energy band from 


16 to 30 k.eV, detected by the same Instrument at 20:55:06 UT. The dominant feature 

IS s 

around footpoint B has an area of appro '•rimately 1.4x10 cm , and emits 30% of the 
total flux. The soft X-ray Image at 20:55:06 UT, simultaneous with the appearance of 
the hard X-ray fuotpolnts, indicates that initially the thermal plasma is denser close to 
the base of the flaring loop system. 

The temporal variation of the turbulent velocity derived from the Calcium spectra 
is compared in Fig. 1 with the hard X-ray emission detected by the SMM Hard X-ray 
Burst Spectrometer (Orwlg et ai. 1980). The emission Is integrated over the energy 
baud from 25 to 386 heV. Significant non-thermal velocities are observed as long as 
the hard X-ray flux remains at high level, and has a hard spectrum, that is, as long as 
there is significant energy release in the flare region. The time variation of the velo- 
city of plasma upflows, as derived from the observed blue shifts, is also shown in Fig 
1 . 

in. A NOM-THERMAL FLAKE LINE BROADENING MODE!,: IViAGNETIC STOCHASTICITY 

How are these data to be interpreted? While there is consensus on the association 
of the excessive line broadening with mas: inoliuns, the nature of these mass motions is 
not as woU ostablishod. In this section, we begin *^y briefly reviewing one common 
explanation for the broadeinog, as set out most recently by Jioschek et al. (1985); our 
aim is to show that this most straightforward model does leave unexplained several cru- 
cial features of the data, which will figure prominently in defining our model. 

a) The Simplest Miniel for Non-thermal line Broadening. 

Let us first consider the most straightforward explanation for the observed broaden- 
ing. As discussed by Doschek et al. (1985), broadened lines may result because of con- 
vective flows within loops, which originate during the "evaporation" phase. This 
Interpretation is suggested by the fact that both non-thermal motions and pleisma 
upflows in the active coronal loops are observed only during the impulsive phase, and 
by the correlation between the peak values of the isotropic non-thermal and upward 
Mow velocities. Furthermore, numerical simulations of flaring loops (see, for example 
Peres et al. 1985) predict line broadening (as well as blue-shifted spectral components), 
even for the case of a single loop configuration; in the case of the single loop, the 
broadening results when the flare heating spatial dependence (Cheng et aJ. 1984) or the 
loop geometry itself is asymmetric, and can also ari.se for certain projections of the loop 
with respect to the llne-of -sight. Moreover, it has been argued that a much more plau- 
sible flare configuration than a single "flare loop" is a complex of loops which 
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participate In the flaring process. In this more realistic case, the line broadening from 
convective flows is yet more easily achieved because the superposition of the contribu- 
tions of a distribution of loops (and their associated convective flows) within a resolu- 
tion element must be taken Into account. In favor of this interpretation is also the fact 
that line broadening decreases with decreasing temperature, which is to be related to 
the lower evaporation velocities expected for cooler material (Doschek et ai. 1985). In 
summary, it seems to us that there is no doubt that Introduction of sufficient geometric 
complexity can allow the above model to account for the data. However, given this 
necessary geometric complexity, it seems worthwhile to ask whether a more ecoiiomlc 
explanation for the phenomena is available. To better define the requirements for such 
an exlanatlon, consider the major observational challenges faced by the geometric 
models: 

1. Lack of a louf^lude dependence of the isotropic turbuient flow parameters. If 
"evaporation" Is the cause of the convective motions which lead to the line broadening 
(as well as to the blue wings), we would expect some dependence of line width on the 
flare longitude (because the flows Induced by the "evaporation" process are dominantly 
vertical to the solar surface). This expected effect would be most pronounced for limb 
f] '-res. It is relatively trivial to .aodel this behavior: Assume tfc ,t the evaporative flow 
within a single loop is described by the blue-shifted component seen disk-on, and 
further assume that the Mare site Is a composite of loops with a uniform distribution in 
the angle between the solar surface normal and the normal to the loop plane (so that 
the loops are on average inclined to the solar surface at angles between 0 and ^ x/2). In 
this most favorable case (because it includes the implausible case of loops lying parallel 
to the solar surface), the convolution of the projected velocity distribution with the 
loop orientation distribution yields a predicted line profile at disk center which will not 
show any red wings (see 2. following), but will show a blue wing. In contrast, unless 
there is a strong asymmetry of flows within loops, a flaring region at the limb uall 
show only weak blue wings {as well as weak red wings). Indeed, the data do show 
that the tr ansien t blue component (which is most naturally associated with the "eva- 
porated" material) is oqI observed in limb flares. If the non-thermal broadening is in 
fact related to these upflows, one would then expect at least some corresponding 
decrease in its contribution to the observed line profile in limb flares (as the systematic 
blue wing is absent): however, the available data do uiH confirm either this effect, or 
the existence of any correlation between the broadening and the angle between the 
line -of -sight and the local radial direction in the immediate vicinity of a flare. Thus, 
the available data do suggest that the turbulent motion field within the ioop(s) is 
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Isotropic, at least to the accuracy which can be attained to date. 

2. Symmetry of the broadened line tor on-disk flares. The dominant changes in 
the Ca XIX line at flare onset are the strong increase in intensity and the substantial 
symmetrical broadening of the line about line center, the latter Independent of flare 
location on the disk. We find it physically implausible to construct an ensemble of 
flaring loops such that the observed Ca XIX line emission receives equal contributions 
from up and down flowing large-scale flows independent of the flare location (as 
required by the line symmetry). Asymmetries in the impulsive heating within any 
given loop can lead to flows of opposite sign near the loop’s two foot points (S. 
Antlochos, private communication), but it is unlikely that such asymmetries always 
arrange themselves so that the emission measure contributions from up and down flow- 
ing material are always identical (to present experimental uncertainties). Antlochos has 
also suggested that during the initial heating phase (before evaporation lias occured), 
the expansion of heated pre-existing coronal material will lead to an initial phase of 
downflows near the loop footpoints (as in fact seen and discussed, for example, in the 
simulations reported by Peres et al. 1985); however, this eflect (which does explain 
the red-shifted portion of the broadening) fails to account fur the comparable blue- 
shifted portion of the broadening. 

8. helative Timing of non-thermal broadening and blue shifts. Non-thermal 
velocities aie often obsei ved before the appearance of high-speed upflows (l.e., before 
the appearance of significant blue wings; Figure 1); certainly, the emission which 
occurs this early during the rise phase cannot be easily explained in terms of a distribu- 
tion of convective velocities driven by the "evaporation" process (see also 2. above). In 
this initial phase, moreover, the non-thermal velocities are quite large; indeed, the asso- 
ciated line broadening (together with a moderate flux of hard X-rays) is often the first 
evidence for flare onset as seen byr looking at high-energy photons. 

In summary, observation *• of individual flares indicate the presence of fairly isotro- 
pic flows during the ear ly part of flare onset; these flows exist In addition to the sys- 
tematic upward flows (which are most likely related to the initial "evaporation" of 
chromo.sphpric material). The lack of any longitude dependence of both the red and 
blue wings of the line broadening for different flares and its symmetry suggests that 
the broadening Is indeed isotropic (and not simply due to a superposition of Doppler- 
shifted components which arise within distinct, but unresolved, loops); the timing of 
the onset of .strong line broadening calls into question its association with the evapora- 
tion process; and the lact that broadening occurs (and is largest) at the very onset of the 
flare argues that It Is pre-vxisUug (not newly-evaporated) coronal material which Is 
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responsible for the line broadeninjs. Otir arguments thus support the proposition that 
the observed broadening appears isotropic not because of accidental surer position of 
various unresolved flaring loops (and associated Interior convective flows) along the 
llne-of -sight, but because the fluid which gives rise to the observed emission is indeed 
Isotropically turbulent in any one given loop structure (l.e., the fluid is, on the dimen- 
sions of the loop, locally isotropically turbulent). 

The idea that the observed broadening is the result of locally-isotropic motions 
within well-defined, local regions in a loop is not new, and in fact underlies the recent 
"maser" model of Melrose and Dulk (1984). They suggest that the flaring coronal 
plasma can be heated to temperatures of ^ 10^ K through local absorption of photons 
produced by direct production of electromagnetic radiation via a gyrosynchrotron loss 
cone instability (viz., Wu and Lee 1979). This instability is thought to occur in a large 
number of small regions, resulting in short elementary bursts; its end result is to 
transform the free energy residing in the perpendicular (to S) motions of mirrored 
(reflected) accelerated electrons to radio frequency radiation. The localized heating in 
the flare region leads then to rapid expansion of the heated plasma at speeds up to the 
ion sound speed; this expansion may account for the observed non-lhermal line broaden- 
ing. Thus, in this model non -thermal motions, and the issociated line broadening, are 
expected as long as the hard x-ray emission is observed. In the next section, we propose 
an alternative picture, which also leads to fairly isotropic motions within any one given 
loop right at flare unsel. 

h) Fisld Lina Stix:hastictty and Tha Onsat ot Flara Fluid Turbulence. 

It has been long recognized that it is very difficult to const uct stable magnetostatic 
(or maguulohydrodynamic) equilibria; that is, it has not in general been possible to dev- 
ise such equilibria (other than for configurations having a high degree of symmetry) in 
a .systematic manner (see., for example, Manheimer and 'ash more- Davies 1984 for a 
review of plasma equilibria and their stability properties). This has led to the conjec- 
ture Uiat this difficulty is not simply computational, but rather is generic, i.e., that 
lack of equilibrium states for geometries that have no simple spatial symmetry is a fun- 
damental property of magnetofluids (Grad 1967; Parker 197P., 1979). This conjecture 
has recently received theoretical support from the demonstration that in the absence of 
rigid boundaries and gravity, symmetric magnetostatic equilibria are topologically 
unstable (Tsinganos 1983; Rosner and Knobloch 1983; Tsinganos, Distler, and Rosnor 
1984); that is, such equilibria have the property that (except ♦'or a special cla.ss of per- 
turbations of zero measure) arbitrary perturbations necessarily lead to a violatrun of the 
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mdguetostatlc constraint - pressure surfaces no longer coincide with majinetlc flux sur- 
faces - so that the system becomes dynamic. Typically, the symmetry of the original 
equilibrium is broken because siugular points arise tArlthln the configuration and, as 
shown by Parker (1982, 1983; see also Vainshtein and Parker 1984), static equilibrium 
is not possible in the vicinity of such new singular points. Thus, the Imposition of 
external perturbations on an equilibrium MHD structure leads to a situation in which 
the system locally departs from static conditions, and local reconnection occurs (cf. 
reviews by Freiberg 1982 and Manhelmer and Lashmore-Davles 1984). In the context 
of laboratory plasma confinement, this process is related to the destruction of magnetic 
surfaces and the formation of magnetic islands; in the context of heating the solar 
corona, Parker (1983) has argued that such a situation arises in the lower solar atmo- 
sphere due to the continual deformation of the coronal magnetic fields by the horizontal 
cellular photospherlc flow field. In summary, local disequilibrium (l.e. lack of local 
static equilibrium) leads naturally to formation of localized reconnection regions and 
Initially Isolated magnetic islands throughout the loop volume. 

What is the consequence of this local destabilization in the context of a single flar- 
ing loop? We conjecture that the system will continue to evolve quasi-steadily (con- 
sistent .vlth the field line topology imnosed by the perturbation) untP island formation 
has proceeded sufficiently to produce regions of overlapping islands; this leads to the 
onset of field line stochastlclty, strongly enhanced local reconnection and dissipation, 
and enhancement of plasma transport coefficients (cf. Freiberg 1982; a related pro<;ess, 
proposed in the context of a specific MHD instability, c.g. Island overlap of unstable 
tearing modes located in distinct unstable tearing layers, has b«3n proposed by Finn 
1975, and applied to the flare onset problem by Spicer 1975, 1976). It is essential to 
distinguish this succession of events from the destablDzlng process previously suggested 
by Low (19B2a,b) and others, who have argued that solar flares represent the terminal 
point of an evolutionary sequence of "nearby" equilibria through which a system 
evolves as it is subjected to external perturbations; this terminal point is reached when 
the sequence of equilibria ends In the sense that there are no longer any nearby equili- 
bria io which the system could evolve. Our model is dlstliiguishcd from these earlier 
studies in that (i) we do not identify the onset of flaring with the point of termination 
of a sequence of nearby equilibrium states, but rather view it as the point at which a 
topological unstable, qi dsl-steadlly evolving state (which need not be force-free; com- 
pare with Heyvaert and Priest 1983) reaches the stochastic domain; and (il) the instabil- 
ity in our picture does not manifest itself in a global MHD instability, but rather in a 
drastic departure from statics throughout the entire volume of the system, wherever 
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singular field line belia\dor as a result of Island overlap, and local reconnection, occurs. 
These locali;£ed reconnection events axe then expected to occur at the very onset of the 
flare -- indeed, in our picture, they mark, the very beginning of flare onset. 

In the same connection, Parker (1983) has followed the hydrodynamics of a layer 

of fluid confin<«l by the pressure of the reconnecting magnetic fields, and found initial 

acceleration of the fluid and formation of Jets that "squirt" away from the reconnection 

site, leading eventually to a quasi-static configuration characterized by jet deceleration 

and simpler magnetic topologies. If this "microscopic" picture of the reconnection 

dynamics applies to the reconnection sites we presume are distributed throughout the 

flare loop volume, we can decide quantitatively whether the observations of line 

broadening are consistent with the our model, and what conclusions regarding the 

nature of the reconnection regions can be drawn. Specifically, consider the May 21 

flare discussed above, for which the necessary data are in hand. Comparison of HXS and 

magnetograph data (the latter obtained at 20:15 UT, courtesy of J. Harvey) shows that 

at the flare site, the mean photospheric magnetic field strength at me of the flare loop 

footpoints (marked A In Figure Id of Hoyng et al.) is * 160 Gauss, covering an area of * 
17 2 

7.4x10 cm . From SXR Images, we estimate that the loop cross-sectional area 

1 ft 

expands to ^ 1.4x'0 cm in the low corona; magnetic flux is conserved, aen a 

reasonable estimate of the coronal field strength is ^ 80 Gauss. From XRP data (taken 

between 20;53-20i55 UT), we estimate an emission measure of * 5x10 cm , so that 

with the (limeusions given above, we obtain an estimate for the luronal plasma density, 

n^ * 3x10^** cm’‘\ Thus, we can finally estimate the Alfven speed in this region; we 
H -1 

find Vj^ 1x10 cm sec . Thi.s compare.s with the observed turbulent broadening 

7 -1 

velocity ot * 2.2x10 cm sec ; as expected, v^ ^ v^. 

Now, if the bruademu^ is indeed due to the superposition of random outflows 

(whoso rate ought to be of the ordor of the Alfven spood) In many reconnection sites, 

then we would expect that ~ v^; this can be true only if the estimate of the plasma 

density in the reconnection site used above Is in error (it is unlikely that our estimate 

of the field .strength is so seriously in error). A moment’s reflection shows why the 

above ostimato for n indeed ought to be in error : we assumed that the emission measure 

is associated with the entire (16 ) coronal volume estimated from the SXH images, 

with a filling factor of unity, even though the model asserts that the emission is in 

fact coming from a large number of small-scale recuunectiou sites scattered throjgliuut 

the volume of the florlug loop. To be consistent, are thus equate v^ and v^, and then 

11 

determine n; we obtain n » H.HxlO cm . This implies that the filling factor is in 
fact ~ 5% of the total volume seen by the SXR in emission. 
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Now, the nows ejected by the localized reconnection regions will take some time to 
dissipate in the low-viscosity corona; in fact, in the absence of tiirbulent mixinf^ 
(which will further contribute to the Isotropization of the turbulent Hows), the shor- 
test relevant time scale for stopping these flows Is the geometrical crossing time: if we 
estimate the lemtih of the entire loop (from footpoint to footi>oint) to be of order 10^^ 
cm, then this time scale is roughly 450 sec. This is consistent with the e-foldiug decay 
time of the observed turbulent velocity v^, r * 400 sec. 

IV. SUMMARY 

In this paper, we have put forth an alternative model for the recent observations 
that solar flare.s are accompanied by line broadening at the very onset of the Hare, i.e., 
earlier than one would aspect any significant "evaporation" cf plasma from the loop 
footpoints. We briefly recapitulate recent theoretical work on the topological stability 
of magnetic structures, which .suggests that departures from static equilibrium should 
be accompanied by the formation of singular points and magnetic islands throughout the 
volume of the plasma, and that when Island overlap occurs, field line stochasticlty and 
localized reconnection result, manifesting themselves by onset of vigorous flows near 
he localized reconnection regi.ns. We suggest that it is th.= superposition of these 
flows associated with local reconnection sites scattered throughout the flare loop which 
leads to the observed line broadening. More particular.'/ we showed that the observa- 
tions aiu quantitatively consistent with a model in whic.i the turbulent broadening of 
high-temperature plasma is entirely due to the outflowing motions from leconnection 
site.s scattered throughout the flaring loop. 

The model we propose to account for the turbulent line broadening has several addi- 
tional attractive features. The must obvious Is that the reccuuectiuu sites we invoke to 
account for the macroscopic plasma streaming will also be sites of localized particle 
acceleration. As dl.scussed In 5 U, hard x-ray bursts are observed precisely during the 
period In which turbulent line broadening is observed; in our model, the fast particles 
which give rise to the bursts may have their origin in the many scattered reconnection 
sites we’ve invoked. In particular, this implies that particle acceleration does not occur 
at just a very few .selected sites within a loop, but rather occurs throughout the loop 
volume. An especially interesting possibility is then that particles accelerated in any 
one such reconnection region continue to be accelurattMl in other reconnection regions 
they encounter as they traverse the loop; such multiple acceleration encounters may 
provide the Kerml process called for by the particle spectra deduced from the SMM hard 
x-ray and gamma ray observations (viz., Ramaty e( el. 1983, and references therein). 
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FIGURE CAPTIONS 


Figure /; Temporal evolution of the impulsive hard X-ray (HXR) emission in the 
energy hand from 25-386 XeV, detected by the Hard X-ray Btirst Spectrometer; of the 
upflow velocity v* and of the turbulent velocity v^ in the soft X-ray (SXR) plasma, 
inferred from the spectra detected in the Calcium channel (3.165-3.231 A) of the Bent 
Crystal Spectrometer. The period of observation of the first spectrum of Figure 2 is also 
Indicated. 

Figure 2: Comparison of two X-ray spectra obtained with the Bent Crystal Spec- 
trometer on SMM before the impulsive hard X-ray increase, at 20:53:02 UT, (a), and at 
the peak of the thermal phase at 21:07:00 UT, (b), of the May 21, 1980 solar flare. 
The smooth curves represent synthetic spectra computed for the given values of electron 
lumperature T^ and Doppler temperature Td* 

Figure 3: Soft X-ray images of the 21®' of May 1980 flaring region obtained in 
the 3.5-8 keV energy band of the Hard X-ray Imaging Spectrometer. The Images are 
deconvolved to aUow for the triangular response of the spectrometer coUimator. The 
crosses represent the positions of the emission peaks measured in the hard X-ray images 
in the 16-30 keV band. The contour levels correspond to 14, 26, 40, 55, 70, 80, 90, 98 
and 1007. of the peak counting rate. 
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High brightness temperature spikes have been observed during a 
radio flare on the M-dwarf flare star AD Leo (Lang et al. 1983). Their 
high brightness temperature (>10^^ K) and circular polarization indicate 
that a coherent radiation mechanism must be responsible for the spike 
emission. The underlying flare emission, which is identified with a lovr 
polarization, gradual component, was found not to be spiky to within the 
200 ms time resolution ul the observations. This note is concerned 
primarily with this non-spiky emission. 

The non-spiky emission is about 100 times more intense than the 
most intense solar flare emission at 1400 MHz, and about four orders of 
magnitude more Intense than a typical solar flare. The brightness 
temperature of the emission is 

T^, - 6.41xl0^‘^(S/f'^)(d/R)^ , 

b 

where S is the radio flux at frequency f in Janskys, f is in GHz, d is 
the distance to the star in parsecs (4.9 pc for AD Leo), and R Is the 
linear size of the source region in solar radii. Taking R to be less 
than the stellar radius (approximately half the radius ot the sun), the 
brightness temperature is found to be Tg > 3xl0‘‘ K. Hence, the 
brightness temperature of the non-spiky emission is also quite high, and 
a nonthermal emission mechanism is required. 


The emission mechanism for the non-spiky emission may be eithev 
incoherent synchrotron radiation or coherent gyrosynchrotron or plasma 
radiation. Independent of whichever mechanism is responsible, the 1400 
MHz observing frequency puts constraints upon the densltv and magnetic 
field strength in the emission region. In order to avoid suppression of 
the radiation by the thermal plasma, the thermal electron density is 
constrained to n < I .24 xlO^^f ^(GHz) ■ 2.4x10^^ cm“^. For synchrotron or 
gyrosynchrotron emission, the electron gyrof requency must be smaller 
than the observation frequency, giving B < 357f(GHz) - 500 G. Also, the 
synchrotron loss time becomes shorter than the 15 min rise time of the 
flare if B > 500 G. Even if the emission mechanism does not directly 


Involve the magnetic field, this condition is required to avoid thermal 
cyclotron absorption cf the radiation. Therefore, the magnetic field 
strength in the emission region is comparable to or not much greater 
than a value that is characteristic of solar coronal loops in active 
regions. 

If the emission is t.ssumed to be incoherent synchrotron radiation, 
the srurce size and magnetic field are constrained by the fact that the 
source size cannot be smaller than the size that is required for it to 
become self-absorbed at the observation frequency. This constrains the 
linear size of the source to (cf. Pacholczyk 1^70, Ramaty and Petrosian 

,/9 1/' S/' 

R 3.3 , 

where is the component of the magnetic field transverse to the 
direction of wave propagation and R, S, f, and d are in the same units 
as above. For the AD Leo flare this gives R S.ABj^' . Hence, for ■ 

100 G, R it 20 stellar radii. Even for as low as 1 G, R ^ 7 stellar 
radii is required. Taking Bj^ * I G, the brightness temperature of the 
emission is Tg 7 >«10^ K and electron energies of at least I MeV are 
required. In summary, for Incoherent emission, either a very large (>10 
stellar radii) "magnetospheric" source region containing mildly 
relativistic electrons and ~100 G magnetic fields is required, or a 
smaller source of low magnetic field strength (Cl G) and fully 
relativistic electrons is required. In e'ther case, the conditions are 
quite distinct from those that characterize the impulsive microwave 
emission from solar flares (cf. Kundu 1965). 

An attractive coherent mechanism for the non-spiky emission, as for 
the spikes, is gyrosynchrotron masering from electrons reflected in a 
magnetic loop (Holman, Eichler, and Kundu 1980). As was noted by Lang 
et al . , second harmonic emission requires a magnetic field strength of 
250 G. An Interesting possibility is that the spiky and the non-spiky 
(which may in fact be spiky on a time scale less than 200 ms) emissions 
are masering at adjacent harmonics, with the non-spiky emission arising 
from the higher harmonic. Alternatively, the radiation may be plasma 
emission at twice the electron plasma frequency or at twice the upper 
hybrid frequency (see Holman 1983 for a review). Emission near the 
plasma frequency Itself is not as likely, since free-free absorption in 
the overlying stellar atmosphere is strong. The magnetic field and 
plasma paramiters that are required by each of these mechanisms are 
comparable to those that are characteristic of flaring regions in the 
solar corona. 

In conclusion, the AD Leo flare emission is clearly not analogous 
to the impulsive phase microwave emission from solar flares. On the 
other hand, the high brightness temperature, presence of superimposed 
spikes, and observing frequency are consistent with the properties of 
solar declmetric Type IV bursts (see Kundu 1965, Kuijpers 1980). The 
coherent emission mechanisms are difficult to distinguish 
observal tonally, since they all produce high brightness tempeiature. 


narrow bandwidth, highly polarized emission. The different mechanisms 
do depend upon different physical parameters, such as magnetic field 
strength for gyrosynchrotron maserlng a electron density for coherent 
emission at twice the electron plasma frequency, and could be 
distinguished li the values of these parameters in the emission region 
could be independently determined. The possibility that the emission 1 
Incoherent synchrotron radiation can be tested with polarization, 
timing, spectral, and high spatial resolution (VLBI) measurements, 
however, since incoherent emission from relativistic electrons would no 
be consistent with a high degree of circular polarization, rapid 
variability would not be consistent with the large source size that is 
required, and narrow-band spectral features would be indicative a 
coherent mechanism. 
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ABSTRACT 

Two 1" X P fields of the Pleiades region, containing 78 cluster members within a limiting magnitude of 14 
mag and centered on two of the most luminous stars of the cluster (20 Tau and 17 Tau) have been observed 
with the Einstein (HEAO 2) Observatory imaging proportional counter. The e.xposure times (~3— 4 x 10^ s) 
and background level give, at the Pleiades distance (- 127 pc), a mean detection threshold of 10** * ergs s“*. 
We have detected one (out of eight) B stars, one (out of 13) A stars, two (out of 10) F stars, 11 (out of 21) G 
stars and six (out of more than 26) K stars. The brightest X-ray source is Hz II 253 (Gl), with L, ss 10*° ’ 
ergs s"*. 

We derive the maximum-likelihood X-ray luminosity functions for the G and K stars in the cluster, and 
show that for the G stars, the Pleiades X-ray luminosity function is significantly brighter than the correspond- 
ing function for Hyades G dwarf stars. The significantly larger number of X-ray bright G stars than K stars 
(even though the Pleiades K stars appear to be relatively rapid rotators), and the lack of detection of M stars, 
suggest that the connection between stellar rotation and the level of X-ray emission is not as straightforward 
as has been herelofoYe thought. 

Subject Iteadinys: clusters: open — stars: rotation — X-rays: sources 


' Also at Harvard-Smithsonian Ccn'.jr for Astrophysics. 


I. INrBODUCTION 

One of the most striking features of stellar X-ray emission 
discovered with the Einstein {HEAO-2) Observatory is the fact 
that stars of the same spectral type and luminosity class emit 
over a wide range. Since the internal structure and the photo- 
sphere of stars of a given spectral type and luminosity class are 
similar, the range of X-ray emission observed can only be 
explained by relating it to other physical stellar parameters 
which do not traditionally figure in placing stars in the H-R 
diagram. 

For example Pallavicini et ai. (1981, 1982), among others (see 
also Waiter 1981), have investigated the dependence of X-ray 
luminosity on stellar rotation rates for a large sample of X-ray 
sources, and have found that most of the variance in the emis- 
sion levels of stars for spectra! types later than F7 can be 
related to the observed dispersion in rotation rates, T1 us they 
find that the X-ray luminosity is quadraiically related to the 
stellar rotation rate or, more sptecifically, that L, 5: (p sin i)^ °, 
where is the X-ray luminosity and v sin i is the equatorial 
rotation speed projected along the line of sight. This result can 
be understood qualitatively by assuming tliat the mechanism 
responsible for heating the outer atmosphere of late-type stars 
to temperatures of several n:illion K is relatco io dissipation of 
energy in magnctic-fucld-dominated regions at the stellar 
surface. The stellar rotation therefore enters in th e above rela- 
tion because in stellar dynamo theory it is one of the important 
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parameters which determine the rate of magnetic-field gener- 
ation in the stellar interior. Similar arguments can be made 
regarding the observed correlation between stellar activity \ 

indicators and Rossby number (cf., Noyes et al. 1984; Schmitt 
et al. 198S; Micela, Sciortino, and Serio 1984; Serio 1984). 

Since stars lose their angular momentum in the course of 
their evolution, it is also reasonable to investigate the depen- 
dence of X-ray luminosity on stellar age. Stern et al. (1981) 
have in fact sui^'cyed the Hyades open cluster, which has an 
age of -^6 X 10* yr, and have found that the median X-ray 
luminosity (LJ of G stars in this cluster is 30 times larger 
than that of the Sun, consistent with a quadratic dependence of 

on stellar age. In contrast, K.u and Chanan (1979). Fcigelson 
and E>e Campli (1981), and Montmerle et al. (1983) have 
observed pre-main-sequence G stars in regions of recent star 
formation, whose median X-ray luminosity (~10*‘ ergs s“‘) 
suggests a linear decay of L, with stellar age; a similar result 
was obtained by Micela et al. (1983) from preliminary analysis 
of X-ray observations of stars in the Pleiades. This paper 
addresses the question of the stellar activity-age relation by 
examining the complete results of an Einstein X-ray survey of 
the central region of the Pleiades, a well-known galactic open 
cluster somewhat younger than the Hyaoes. 

Our paper is organized as follows: we describe the data 
selection procedures, the method of analysis, and the results of 
the survey in the next section. The implications of these results 
for the dependence of X-ray emission from late type stars on 
the properties of the underlying star are discussed in § III. Our 
results are summarized in § IV. 

II. DATA SELECTION, METHOD OF ANALYSIS, AND RESULTS 

Two Einstein IPC (Gorenstein, Hamden, and Fabricant 
1981; Giacconi et al. 1979) exposures centered on two of the 
brightest stars of the Pleiades (17 Tau and 20 Tau) were 
obtained on 198! February 7 and 8. The two fields overlap by 
approximately 50% in the ~60' x 60' useful region of the IPC 
field of view, and consequently 40 Pleiades stars fell into the 
overlap region. The characteristics cf the two exposures are 
summarized in Table 1, while Table 2 lists those stars which are 
included in the field of view of the two exposures which are 
part of the cluster according to the Hertzsprung catalog (1947, 
hereafter Hz 11). At the distance of the Pleiades, this catalogs 
completeness limit of apparent magnitude 14 correspond to 
spectral ty late K. 


a) Analysis Procedures 

The data have been analyzed by means of the revised Ein- 
stein IPC software (Hamden et al. 1984).* Images were search- 
ed for 3 a fluctuations with respect to the local background in 
the energy band 0.2-3.5 keV. Vignetting, detection-cell exten- 
sion, and shadowing by the IPC entrance aperture ribs were 
taken into account to give the effective count rate for each 
detected source together with its statistical error. With the 
detection threshold at 3 a, approximately two of the sources 
are expected to be spurious. The positions of the detected 
X-ray sources were compared with the positions of all the stars 
(with m,. < 14) falling in the useful field of view. We have used 
as identification criterion the condition that the position of the 
X-*ay source be within 1' of the position of the optical star. 
This procedure gives a probability of chance identification with 
cluster members of lc‘j than 5 x 10'*. 

In addition i searching for X-ray sources in the field, we 
have also obtained count rates (corrected for instrumental 
effects as above) at the locations of all cluster stars in our 
catalog. Table 3 summarizes our results: We provide lumin- 
osities for all stars detected above the 3 a level; call out all stars 
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for which the detection significance level lies between 2 a and 3 

a, and give the corresponding luminosities; and provide 2 a 

upper bounds for the remaining stars. Random fluctuations are ^ ^ ^ 

expected to occur in -^2% of the trials in composing such a 

list; hence, about two false sources appear in the 2 a portion of 

Table 3. for a total of perhaps four false sources in Table 3.^ 

^ One of the imponant differences between this retised software and the 
former version is in the source-detection algorithm, which searches for fluctua- 
tions with respect to the local background rather than with respect to an 
average background. This ensures more reliable results for sources near the 
detection threshold, as well as in crowded fields. 

’ The results of a survey of the same region reported by Caillault and 
Helfand 1985 are not directly comparable to ours for the marginally detected 
sources, since they used the earlier “ Rev O ” Einstein software. 

Wc have converted from count rates to X-ray luminosities 
using a conversion factor of 2.2 x 10"“ ergs count'* cm'*. 

This value was obtained by assuming a source temperature of 
~3 X 10*’ K and a hydrogen column density of 2.0 x 10^® 
cm'^; the latter was computed according to Bohlin, Savage, 
and Drake (1978) using a reddening of £(B— K)%0.04 
(Crawford and Perry (1976). The total conversion factor for 
count rate to luminosity, assuming a common distance of 1 27 
pc for all cluster members, is then 4.35 x 10^' ergs count'*. 

h) Results 

We show a map of the region surveyed, giving the locations 
of both the Pleiades stars in the two fields and the detected 
sources, in Figure 1. The inferred X-ray luminosity along with 
its statistical uncertainty, or the 2 a upper limit, are listed in 
Table 3 for all stars in our catalog, as explained above. The 
errors listed in Table 3 are the statistical errors. The estimated 
overall error in the X-ray luminosities is ~40%-60%; if con- 
sists of statistical errors (ranging from 15% to .''•'%), systematic 
errors in instrument calibrations (<20"o; Hamden ef al. 1979), 
errors in the individual cluster member distance (<3%), and 
systematic errors in converting counts to flux due to the 
assumed hydrogen column density and source temperature 
(< 20 %.)- 

We have detected a total of 21 distinct X-ray sources, associ- 
ated with cluster members. 14 of which are significant at above 
the 3 a level and 7 at above the 2 a le\el. More specifically we 
have detected: one (1 -f 0) out of eight B stars; one (I 4- 0) out 
of 13 A stars; two(l -i- 1) out of 10 F stars, 11(11 -t- 0) out of 21 
G stars; and six (1 -I- 5) out of 26 K stars in the Hz catalog 
(since the Hz catalog is not complete for apparent magnitudes 
greater than 14, this corresponds to an upper limit on the 
fraction of K stars detected). 

Three additional X-ray sources detected at the 3 a level do 
not correspond to any catalogued Pleiades star within a dis- 
tance of 1' from the detection-cell center. Two of these X-ray 
sources have been identified with field stars using the same 
identification criterion as discussed above; the remaining star 
may well correspond to an uncatalogued late K or M star. 

Our results are summarized in graphical form in Figures 
2-4. Figure 2 shows an H-R diagram of the cluster stars in the 
two fields of view for which count rates were obtained. The 
statistics of the breakdown into spectral types of the detected 
stars is given in Figure 3; Figure 4 summarizes the inferred 
X-ray luminosities and upper limits obtained in this survey. As 
can be seen from Figure X we have values foi the X-ray lumin- 
osity for more than 50% cf the G stars in our sample; hence 
the median of the X-ray luminosity function for G stars in the 
Pleiades can be estimated from Figure 4 as the luminosity of 
the weakest X-ray source identified with a G star: 

(median) = 10^® * ergs s'*. 
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'll. DISCUSSION 

In this section we discuss the conclusions which can be 
drawn from our survey, with particular attention to the depen- 
dence of the level of stellar X-ray emission on age. 

a) Comparison with Other Surveys 

It is particularly useful to compare the data presented here 
with the results of the Einstein Hyades survey (Stern et al. 
1981). Only six stars (~7% of the stars from A to K) in the 
Hyades survey have been detected as X-ray sources with 
luminosities above the threshold for our Pleiades survey. Since 
the primary difference between these two open clusters is their 
age (the Pleiades being 1/10 the age of the Hyades), it is 
natural to attribute this difference in X-ray luminosity to this 
age difference (or to stellar attributes which vary with stellar 
age); this point of view is also consistent with the expected 
decline of stellar activity with time for any given star. 

More specifically. Stern et of. (1981) detected ~80% of the 
Hyades G dwarf stars, finding that the median X-ray lumin- 
osity for these stars (-^10^® ‘ ergs s"’) is approximately 30 
times the X-ray luminosity of the Sun (which is 10 times 
older than the Hyades stars). Among the Hyades stars, three 
out of 13 G dwarfs have X-ray luminosity above of 2.5 x 10^® 
ergs s"'; in comparison, we find that more than half of the 
Pleiades G dwarfs emit above this threshold. Since the Hz 
catagog is not complete for K stars, the comparison of detec- 
tion statistics for K stars is not as meaningful as for the G stars. 
However, while Stem et al. did detect in X-rays only one out of 
the 18 K stars in their combined field of view at a level above 
3 X 10^® ergs s"‘, we have detected six (or seven, if one 
assumes the source which has not been identified to be a K 
star) out of more than 26 K stars. Thus a clear decline in 
activity level (as measured by the stellar X-ray luminosity) with 
age is evident. This decline is illustrated in Figure 5, in which 
we plot the average detection threshold, the median, and the 
highest detection for X-ray observations of G spectral type T 
Tauri stars (Ku and Chanan 1979), Pleiades stars (this paper), 
Hyades stars (Stern et al. 1981), and local disk population 
members (N'aiana et al. 1981). 

In order to quantify this behavior we have constructed 
maximum-likelihood luminosity functions (Avni et al. 1980; 
Schmitt 1984) for both the Pleiades and Hyades G dwarf stars. 
The Pleiades data used for these calculations are given in Table 
3; the rc.sults of our computation for the Pleiades and Hyades 
are shown in Figure 6. We have applied a “bootstrap” calcu- 
lation in order to compute the mean luminosities for the two 
clusters, including errors (Egron 1982; Schmitt 1984); we find 
that 

= 10” ' ergs s-‘ ; 

cigss-*.-* 

4 We noie. however, lhal ihe results reported by Stern ei al. (1981) were 
obtained with the "Rev O" version of the Einswin data-analysis software, 
w'hile our results for the Pleiades stellar cluster we obtained using the updated 
“ Rev I " software; for this reason, there might be small ( ~ 10%) discrepancies 
between our paper and theirs in the luminosity levels of stars. 

We can therefore conclude that the Pleiade.« G stars have a 
significantly (>3 a) higher level of X-ray emission than the 
Hyades G stars. 

We have also compared the range of the X-ray luminosity 
for Pleiades stars in our survey that have 0.4 < fl — F < 1.0 
with a similar sample of stars belonging to the Ursa Major 
cluster (Walter et al. 1984); specifically, the value of L, for 
UMa members ranges between 10^" ® and 10^®® ergs s“‘, 
while for Pleiades stars in our survey L„ range; between 10^® ^ 
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and 10^° ^ ergs s~‘. Because the UMa cluster is older than the 
Pleiades cluster, this result is consistent with a trend of decreas- 
ing level of X-ray emission with stellar age. A more detailed 
comparison is not possible at this time because of the lack of 
upper bounds information for the UMa cluster. 

In spite of the above reassuring evidence that stellar activity 
indeed declines as the stellar population ages, the Pleiades data 
also suggest that the sample picture in which the stellar- 
activity level of late-type stars is dominantly determined by the 
stellar-rotation rate, so that the age dependence is introduced 
by, for example, slow-down in the stellar spin rate, is not the 
entire story. In particular the absence of sources definitively 
identified with M dwarf stars is very puzzling. To see this, note 
that the median X-ray luminosity of Pleiades dG stars is SO 
times as large as the value found by Vaiana et a/. (1981) for field 
dG stars. Then, if we assume that the X-ray luminosity function 
for the Pleiades dM stars is shifted by the same amount from 
that of the field stars, assume that the number of Pleiades M 
stars is ^ 10 times larger than that of G stars, (as is roughly the 
case for galactic-disc stellar-space densities; Allen 1976), and 
use the X-ray luminosity function of field dM stars derived by 
Rosner ei al. (1981). we find that ~ 10 dM stars should have 
been seen above the threshold of our survey. Instead, we have 
observed only or.i source (X-3 in Table 3) which could possibly 
correspond to an uncatalogucd M star; there are no other 
unidentified sources. We therefore conclude that the lumin- 
osity scaling with age must be spectral-type-depcndcnt. as 
would occur if there were another stellar parameter (dependent 
upon spectral type) which figures in determining stellar activity 
level. The most plausible such parameter is the outer stellar 
convection depth; for example, Schmitt et al. ('985) have sug- 
gested that the activity level of dF stars is correlated with the 
Rossby number at the basis of the convection zone, and is thus 
related to a particular combination of rotation rate and con- 
vection zone depth (see also Noyes et al. 1984 for such evidence 
derived from Ca il data). We further note that the above argu- 
ment derived from the (non)observation of dM stars is reen- 
forced by a comparison of the X-ray luminosity functions of 
Pleiades dG and dK stars: using the nonparametric tests 
described by Schmitt (1984), we find that they are different at a 
significance level higher than 98%. with the dG star X-ray- 
luminosity function cxiendiFig to higher values of L, than that 
of the dK stars. This difference is quite remarkable in its own 
right and will be discussed further below. 

h) X-ray Luminosity anil Rotation 

Only eight of the stars later than F7 (specifically: four K. one 
G, three F) in our sample have known values of projected 
equatorial velocities i' sin i; for 1 1 K stars, we have upper limits 
to r sin i and for two stars (one K, one G) we have lower limits 
to V sin I. In particular, two K, two G, and one F stars detected 
as X-ray sources have projected equatorial velocity v sin f in 
excess of 45 km s" ', while two K stars with r sin i > 50 km s" ' 
were not detected (although their X-ray luminosity upper 
limits are quite high and hence not very constraining). 

Figure 7 shows a log-log plot of X-ray luminosity and i; sin i 
for stars of types F7 to K in our sample for which velocity 
measurements or upper/lower limits are available. As a com- 
parison we have plotted the best fit for field stars obtained by 
Pallavicini et al. (1982); we note that the stars in the Pleiades 
have X-ray luminosities up to 2 orders of magnitude lower than 
RS CVn-tyjDe stars with comparable rotational velocity, pi 
than the prediction of the Pallavicini et al. (1982) best fit. The 
RS CVn’s are likely to have substantially different behavior in 
activity-rotation correlations from that of the present sample 
of stars because the former are close binaries; therefore we 
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have applied nonparametric maximum-likelihood linear reg- 
ression analysis (including the upper bounds; Schmitt 1984) to 
our data shown in Figure 7 and to the data of Pallavicini et ^ 

(1982), but excluding the RS CVn stars. We obtain log L, % ^ 

27.5 -f 1.5 log r sin i (solid line in Fig. 7), instead of the quadra- 
tic relationship between X-ray luminosity and projected equa- 
torial velocity obtained by Pallavicini et ai (1982). Our 
regression function is, however, a poor fit to these data, sug- 
gesting that a single power law of the type applied here is not 
an adequate description of the data even in the RS CVn- 
excluded combined sample. 

c) The Rapidly-rotating but Relatii^ly Inactive K Dwarfs 

Given the currently popular notion that stellar activity levels 
ought to correlate well with stel!,'>’’ rotation rates, it seems 
surprising that the Pleiades K dwarfs have systematically lower 
activity levels than the corresponding G dwarfs (shown most 
succinctly by the detection histogram in Fig. 3), because the dK 
stars in this cluster include a number of rapid rotators (cf. 

Stauffer et al. 1984). We believe that one can, however, account 
for this behavior. As suggested by the sciar rotation evolution 
models of Endal and Sofia (1981), rapid evolution on the radi- 
ative track leads to spin-up of the entire star, which is followed 
by rapid spin-down (by magnetic braking) of the outer convec- 
tive envelope. If, on average, the K stars in our sample corre- 
spond to the former evolutionary stage, while the G stars Ogain 
on average) correspond to the latter stage, then one would 
expect the rapidly rotating K dwarfs to have relatively modest 
rates of internal differential rc ation at the convection zone- 
radiative core interface, whereas more slowly rotating G 
dwarfs would still have (at least just after the rapid envelope 
spin-down period) a rapidly spinning core and hence substan- 
tial differential rotation at this interface. 

This scenario is supported by the fact that, at the nominal 
Pleiades age of ~ 10^ ’ yr, a solar-like star would have just 
reached the zero-age main sequence, as shown by Figure 5 of 
Endal and Sofia (1981). Such a star would already show a 
greater than threefold difference between its surface rotation 
rate and the (larger) radiative core rotation rate. In contrast, a 
lower mass K dwarf at the Pleiades age would not have 
reached this stage as yet. and the consequent maximum differ- 
ence in rotation rate between the core and the surface would be 
less than twofold. 

As has been recently suggested by Rosner(1980), Spiegel and 
Wei't (1980), and Golub et al. (1981), the presence of a shear 
near the base of an outer convective zone is highly favorable to 
a "shell dy.”imo'’ driven by a gradient in s2 at the convective 
envelope-rauiative core boundary. Hence, one would expect 
on this model, at least qualitatively, that magnetic activity is 
most vigorous just at the point at w hich this differential rota- 
tion is maximized, thereby accounting for the observed differ- 
ence in activity levels in the Pleiades G and K dwarfs. This 
scenario would also account for the paucity of X-ray sources 
associated with the dM stars in our Pleiades fields. Tbis inter- 
pretation is not in contradiction with the observed high level of 
X-ray emission in T Tauri stars (e.g., Ku and Chanan 1979) 
both because the interior state of these stars is highly uncertain 
and because of the possible influence of primordial ma<'nctic 
fields. 


d) A Stars 

The only A star which has been detected in X-rays in this 
survey (Hz II 1384) has an extremely high X-ray luminosity: 
L, 5: 1.3 X 10*° ergs s"‘, i.e.. three orders of magnitude higher 
than the value indicated as typical by Golub et al. (1983) in 
their survey of A stars detected by Einstein. As shown by 
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Schmitt et al. (1984), A stars detected at peculiarly high values 
of X-ray luminosity are very likely in multiple systems in which 
the X-rays are dominantly emitted from an unseen late-type 
companion; however, Abt et al. (1965) have not found Hz II 
1384 to be a spectroscopic binary. Furthermore, if we are to 
take the argument presented above seriously (and consid'^r the 
relative dearth of detections of such stars), then any dK o: later 
spectral-type companion ought not to be a vigorous X-ray 
source; instead, the X-ray luminosity of Hz II 1384 places it at 
the very upper end of the dG star luminosity function, and 
hence will above the emission levels observed and expected for 
later spectral type stars. Thus, although there is substantial 
observational motivation based on surveys of held dwarf A 
stars to place a strong upper bound on the possible contribu- 
tion of the A star (in the hypothesized binary system) to the 
observed X*ray luminosity, the alternative (namely the pre- 
sence of any lower-mass companion to this dwarf which is 
responsible for the observed emission) appears to be problem- 
atic (as It would require a dG companion which is exception- 
ally bright in X-ray and which has escaped optical detection). 

Given this difficulty, we would like to consider one remain- 
ing alternative explanation, which is even more speculative. A 
stars also go through a fully convective (well-mixed) stage 
during the descent to the main sequence; in this well-mixed 
regime, large radial gradients in rotation rate are unlikely. 
However, since the radial contraction which occurs during 
descent to the main sequence is nonhomologous, there must be 
a regime during which the star is (at least partially) convective 
and strongly differentially rotating. These circumstances are 
suitable for magnetic dynamo action, but continued dynamo 
action is not very plausible once the star has reached the main 
sequence (because convection ought to have ceased). Hence it 
may be that the observed activity is due tc the dissipation of 
the magnetic fields produced during such an evanescent 
niagnetically-active regime. This speculative suggestion seems 
to be consistent with the fact that the normal-field d.4 stars 
show no evidence for “activity” (Schmitt et al. 1985), but that 
dA stars known to be associated with strong surface magnetic 
fields (i.e.. .4, stars such as tv Oph and a CVn) have been 
detected in X-rays, with L, =e 10*’ ergs s" ' (Golub et al. 1983). 

e) I'uriabiliiy 

We have searched for variability between the two IPC obser- 
vations (i.e.. on a time scale of ^ 1 day) for all the stars in the 
overlap region which were detected above the 3 a level in both 
exposures. Our analysis detected variability (at a level of 3.3 o) 
only in the X-ray flux of the star Hz II 303, which varied by a 
factor of ~ 3. No significant variability information could be 
obtained for the remaining stars in the overlap region. 

IV. SUMMARY 

We have analyzed the data obtained with two pointed obser- 
vations of 1“ X 1‘ fields of the Pleiades region, have derived the 
maximum-likelihood X-ray luminosity functions for the 
Pleiades G and K stars in the cluster, and have shown that for 
the G stars, the Pleiades X-ray luminosity function is signifi- 
cantly brighter than the corresponding function for Hy.j(tes G 
dwarf. This finding indicates a dependence of X-ray luminosity 
on stellar age, which is confirmed by comparison of the same 
data with median X-ray luminosities of pre-maiu-sequence 
and local-disc-population dwarf G stars. Furthermore, 
although they arc bright in X-rays, the X-ray luminosity of 
late-type Pleiades stars falls generally below the predictions of 
(he L, it (v sin i)^ relation of Pallavicini et al. (1981), for late- 
type main-sequence stars, indicating that such a simple relation 
between X-ray luminosity and rotational velocity for late-type 
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stars is not an adequate description of the full data. 

We have suggested that the significantly larger number of “7 

bright X-ray sources associated with G stars than with K stars, 
the lack of detections of M stars, and the relatively rapid rota- 
tion of the Pleiades K stars (Stauffer et at. 1984), can be 
explained in terms of the onset of internal differential rotation 
near the convective envelope-radiative core interface afte the 
spin-up phase during evolution to the main sequence. In this 
picture, the K stars are still spun-up, and have relatively little 
differential rotation, while the G stars have already started the 
spin-down by magnetic braking, thereby developing a gradient 
in angular velocity at the boundary between a fast rotating 
core and a slower convective envelope. Finally the detection of 
an A star in the cluster as a vigorous X-ray source cannot be 
readily accounted for by invoking either intrinsic emission 
along the lines of the '‘solar analogy," or an unresolved low- 
mass. active companion, raising the question whether the 
observed emission is the “aftermath” of an earlier epoch of 
magnetic dynamo activity during the descent of the 
(convecting) progenitor star to the main sequence. 
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TABLE I 

CiuRAf rtRisTirs or th£ Ob»rvation5 


I-S457 1-5458 


Center coordinain (1.^1 3*4I"'54M .1*’42"'5(y8 

23 57 27:8 24 12 46:9 

Start lime (JD) 2142.8131 2143.9195 

Effective exposure (SI 2931.0 4734 4 

Gusier stars (m, < 14) 57 61 

Limiting sensitivity ~2.8 x 10'* ergs s"' -2.3 x 10'* ergs s'* 





TABLb 2 

PuiAUU Stars in this Surves' 


Hz II 

R.A (1950) 

Dec! (1950) 

m, 

Ref 

B-V 

Sp.* 

Ref. 


Ref. 

P‘ 

Notes 

120 

3‘40"34M 

23 3103" 

10.79 

1 

070 

F6 






158‘ 

3 40 44.5 

24 1306 

3.23 

1 

0.25 

A7 

2 

70 

3 



193" 

3 40 52.2 

24 05 29 

11.29 

1 

0.81 

G7 


. . . 


097 


232 

3 41 01.3 

24 24 04 

8.06 

1 

0.20 

AS 

2 

20 

3 


HD 23194 

233 

3 41 00.7 

23 43 37 

9.66 

1 

0.52 

n 


80 

3 



253 

3 41 04.6 

24 20 54 

10.66 

1 

0.68 

Gl 

4 



. . 


263* 

3 41 06.2 

24 07 10 

11.54 

1 

0.88 

G8 

4 



0.98 


298* 

341 144 

23 52 34 

10.86 

1 

0.88 

G 4- K 




.1. 

e 

299* 

3 41 14.7 

23 52 30 










mJ3* 

3 41 16.2 

23 5647 

10.48 

1 

0.89 

G9 






314 

3 4 1 20.8 

24 38 26 

10.56 

I 

064 

Gl 

4 





320 

3 41 21.2 

24 3 7 02 

1104 

1 

0.88 

G5 

4 

. . 




324 

3 41 22 6 

24 36 46 

13.00 

5 

1.07 

K2 


90 


... 


335* 

3 41 24.6 

23 54 46 

13.76 

5 

1 28 

K5 

6 



0.56 

f 

338 

3 41 25.1 

23 58 37 

9.07 

1 

046 

F2 

2 

<40 

3 

. . 


344* 

3 41 26.9 

24 14 21 

8.17 

1 

0.27 

A8 

2 

200 

3 

. . 


345* 

3 41 27.2 

24 2603 

11.65 

5 

0.85 

G8 

4 

. . 


0.99 


357 

3 41 29 4 

24 00 58 

13.32 

7 

1.19 

K6 

6 

<10 


0.99 

f 

405 

3 41 41.5 

24 39 47 

9.83 

1 

054 

F9 

2 

15 

8 

. . 


430* 

3 41 45.3 

24 04 33 

1145 

5 

0.85 

G8 

4 



0.99 


447* 

3 41 49.5 

24 08 04 

546 

I 

-0.04 

B7 IV 

2 

260 

3 


16 Tau 

468* 

3 41 54.0 

23 57 29 

3.71 

1 

-Oil 

B6 III 

2 

220 

3 


17 Tau 

476* 

3 41 55 6 

23 45 58 

1081 

1 

0.80 

Gl 

4 



... 

HD 23302 

489* 

3 41 57 6 

24 16 39 

10.36 

1 

063 

GO 

*> 4 





522* 

3 42 05.1 

23 41 05 

1197 

5 

092 

K2 

9 

. . . 


0.99 


530 

3 42 07.2 

23 32 52 

8.95 

I 

0.39 

F3 


<40 

3 

. . 


531* 

3 42 07.9 

24 06 31 

858 

1 

0 31 

Am? 

2 

75 

3 

. . 


541 

3 42 10.4 

24 41 03 

5.65 

1 

-0.07 

B8 


245 

3 

. . 

HD 23324 

554 

3 42 12 9 

24 25 53 

14.09 

9 

092 

K5 

9 



0.94 


563 

3 42 13.5 

24 18 45 

4.31 

I 

-Oil 

B6 

2 

135 

3 

... 

19 Tau 

625 

3 42 23.1 

23 34 23 

12.57 

I 

1 10 

K3 


>50 


096 

8 

652 

3 42 27.7 

23 52 50 

8.04 

I 

021 

A3 

2 

235 

3 



676 

3 42 31.4 

23 36 21 

13.71 

5 

1.31 

K3.5 

6 

<10 


O.W 

f 

686* 

3 42 34.2 

24 08 56 

1362 

1 

1.04 

K2 


150 


099 

fg 

697* 

3 42 35.5 

24 18 32 

860 

1 

035 

A9 

2 

75 

3 

... 


708* 

3 4 2 36 9 

23 5544 

1013 

1 

062 

GO 

4 

70 

8 



717* 

3 42 39.0 

24 1 0 53 

7 18 

1 

016 

AI 

2 

15 

3 


HD 23387 

727 

3 42 41 1 

24 28 23 

9’0 

1 

0 55 

F9 

•> 

45 

8 



738 

3 42 41 3 

23 35 60 

12.26 

1 

1 16 

G5 

to 

>60 


0.93 


745* 

3 42 42 7 

24 08 03 

9.45 

1 

0.52 

F5 

2 

65 

8 

... 


746* 

3 42 42.9 

24 16 38 

11 27 

1 

0 92 

G5 

4 

... 


0.99 


761* 

3 42 45.7 

24 03 58 

10.55 

1 

0 67 

Gl 

4 





785* 

3 4 2 50 8 

24 12 49 

3.88 

1 

-0.07 

B7 III 

2 

40 

? 


20 Tau 

7V3* 

3 42 50.7 

23 41 56 

14.305 

5 

I 40 

K8 


<10 


0.96 

f 

’99* 

3 42 52.2 

23 43 12 

13.71 

5 

1 31 

K3 

10 

<10 


0.96 


804 

3 42 53.2 

23 5305 

7 85 

1 

0.20 

A2 

2 

170 

3 



817 

3 42 55 4 

24 2402 

5.76 

1 

-004 

H8 

2 

220 

3 

... 

21 Tau 

859 

3 43 03.9 

24 22 26 

6.43 

1 

-002 

B9 

2 

250 

3 


22 Tau 

870 

3 43 04 7 

23 3500 

12.72 

1 

1.24 

K5 




0.93 


879* 

3 43 07.4 

24 24 49 

12.82 

5 

108 

K2 


<10 


0.96 

e. g 

883* 

3 43 07.8 

24 24 33 

13.05 

5 

1.15 

K4 




0.99 

e. f 







TABLE 2 — Coniinutd 


Hz II 

R.A. (1950) 

Decl (1950) 

"•i 

Ref. 

B-y 

Sp* 

Ref 


Ref. 

r 

Notes 

885 

3'43*08‘3 

24'42'47* 

1105 

1 

1.01 

K3 

4 

<10 


0.97 


916' 

3 43 12.6 

24 28 07 

11.71 

1 

0.87 

G8 

4 



0.99 


956' 

3 43 17.3 

24 02 10 

7.96 

1 

0.32 

ro 


150 

3 

.. . 

HD 23479 

974 

3 43 21.1 

24 37 55 

13.86 

5 

1.32 

K7 

9 



0.95 


980 

3 43 21.1 

23 47 41 

4.18 

1 

-0.06 

B6 IV 

2 

2'75 

3 


23Tau 

996' 

3 43 23.5 

24 2500 

1042 

1 

0.64 

Gl 

4 



... 


1028' 

3 43 28.5 

24 0606 

7.35 

1 

0.10 

A2 

2 

110 

3 

. . 


1032' 

3 43 29.4 

24 16 50 

11 34 

1 

0.86 

G8 

4 

... 


0.99 


1061' 

3 43 32.5 

23 57 49 

14.28 

5 

1.42 

K5 

10 

<10 


0.95 

f 

1094' 

3 43 37.5 

23 48 49 

14.02 

1 

140 

K8 




0.91 

e, f 

1100' 

3 43 38.4 

24 11 24 

1116 

1 

1.15 

K3 

4 

<10 


0.92 

f 

1110' 

3 43 39 8 

24 2201 

13.41 

9 

1.23 

K6.5 

6 



0.99 


1117 

3 43 39 5 

23 38 

10.20 

1 

0.73 

G6 

2 





1122 

3 43 40 7 

23 56 60 

9.29 

1 

0.46 

F4 

2 

28 

10 



1124' 

3 43 40.8 

23 52 35 

1112 

1 

0.92 

KO 

4 

<10 


0.97 

8 

1207 

3 43 55 4 

24 38 36 

1047 

1 

062 

Gl 

4 

... 


... 


1266 

3 44 04.1 

24 4001 

8.28 

1 

0 36 

F2 


95 

3 



1280' 

3 44 04.9 

24 00 25 

14.55 

5 

1.37 

K7 

10 



0.97 


1284' 

3 44 05.7 

23 50 32 

8 37 

1 

030 

A9 

2 

100 

3 

. . 


1298 

3 44 08.6 

23 33 44 

1118 

1 

1.02 

K2 


10 


0.98 


1355 

3 44 19.6 

23 5303 

14 07 

7 

1.40 

K5 

6 

<10 


0.77 


1380 

3 44 22.6 

23 3902 

699 

1 

003 

Al 

2 

235 

3 

... 

e 

1384 

3 44 24.8 

24 26Ci9 

7.66 

1 

021 

A4 

2 

215 

3 

. . 


1431 

3 44 30 6 

24 08 09 

681 

1 

006 

A40 

2 

40 

3 

... 


1454 

3 44 34 3 

24 31 54 

12.82 

1 

1 14 

K5 

10 

<10 


0.97 


1514 

3 44 34 4 

24 12 44 

10.48 

1 

064 

Gl 




... 


1516 

3 44 4) 4 

24 08 58 

13.87 

5 

1.27 

K6 




0.99 

f 

1531 

3 44 42 9 

23 49 1 1 

1341 

8 

1.22 

K5 


50 


0.94 

g 


‘ Sptciral lypt computed according to Allen 1976 from the quoted B-\' values, unless otherwise noted. 

* Equatorial velocti) projected along the line of sight (km s"'(. from Stauffer or al. 1984. except when otherwise noted. 

* Probability of duster's membership according to Jc lies 1973. 

* Star present in the overlap region of the two exposures. 

‘ Double system, data from Binnendiik 1946. except when otherwise noted. 

' Hare star; data from Binnendijk 1946. except when otherwise noted. 

' Photometric periodic variable (Van Lcewen and .Alphenaar 1982. Alphcnaar and Van Leewen 1981). 

REnRESCES. — (II Johnson and Mitchell 1958 (2) Mendoza 1956. (3) Anderson. Stoeckly, and Kraft 1966 ( 4) Wilson 1963 (5) Landolt 1979. 
|6lKrafi and Greenstein 1969. (7| Stauffer 1980. Stauffer el uf 1984. (8) Kraft 1967 (9) Jones 1973.(10) Herbig 1962. 
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table 3 


X-Ray Ll'minosihes or Stars in the Survey 


no” ergss-') 

Hz 11 

Spectral type 

I 5458 

I 5457 

120 

F6 


<7.8 

158 

A7 

<2.7 

<4 6 

193 

G7 

6.3 ± 1.6 

3.5 ± 1.2 

232 

A5 

. . . 

<3.2 

233 

F7 


<2.0 

253 

G1 


22.4 ± 3.2 

263 

G8 

4.6 ± 17* 

54 ± 16 

298 

G -E K 

2.9 + 1 4* 

4.8 ± 1.4 

299 

G -E K 



303 

G9 

14 4 ± 2.3 

5.2 ± 1.6 

314 

G1 

<8.7 


320 

G5 

8.8 ± 2.6 


324 

K2 

<8.7 


335 

K5 

<6.i 

<2.6 

338 

F2 


<2 9 

344 

AS 

<1.7 

<2.8 

345 

G8 

12.7+ 1.9 

10.9 + 2.7 

357 

K6 

. . . 

<3.0 

405 

F9 

<44 


430 

G8 

<16 

<2.1 

447 

B7 IV 

<1.4 

<1.9 

468 

B6 111 

<1.7 

<1.5 

476 

G1 

<2.7 

<2.1 

489 

GO 

<18 

<5.2 

522 

K2 

<2.2 

<2.5 

530 

F3 


<2 6 

531 

Am? 

< 1.4 

<1.9 

541 

BS 

<2.9 

... 

554 

K5 

2 0 + 0.9* 

<3.1 

563 

B6 

<2.1 


625 

K3 


59+17 

652 

A3 


<19 

676 

K3 5 


<3.6 

686 

K2 

<20 

3 1 ± 15* 

697 

A 9 

V 1 b 

<3 1 

708 

GO 

47+1.2 

2.9+ 1.2- 

717 

A! 

<1.6 

<2.0 

727 

F9 

9.1 ± 1.5 


7T8 

G5 


6.5 ± 2.0 

745 

F5 

<19 

3.4 ± 1 5* 

746 

G5 

<14 

<3.1 

761 

G1 

3.3 ± 1.0 

<2.7 

785 

B7 111 

<1.6 

<2.4 

793 

K8 

<3 8 

<6 2 

799 

K3 

2.9 ± 1.5* 

<3.7 

804 

A2 


<1.9 

817 

B8 

<1.8 

.. . 

859 

B9 

<16 

... 

870 

K5 


<2.5 

879 

K2 

<18 

<6.0 

883 

K4 

<1.8 

<6.4 



TABLE 3 — Continued 


•Xhl 


Hz II 

Spectral type 

1 5458 

1 5457 

885 

K3 

<2.4 


916 

G8 

<1.5 

<3.9 

956 

FO 

<4.1 

<8.3 

974 

K7 

<2.2 


980 

B6 IV 

4.2 ± 1.6* 


996 

GI 

<2.1 

<4.8 

1028 

A2 

<1.4 

<2.4 

1032 

G8 

8.2 ± 1.4 

6 9 ±2 4* 

1061 

K5 

<2.2 

<30 

1094 

K8 

<2.8 

<3.5 

1100 

K3 

2.8 I 1.0* 

3.7+ 1.7* 

1110 

K6.5 

< 1.6 

< 5.4 

1117 

G6 

. . 

<12.8 

1122 

F4 


<3 5 

1124 

KO 

<3 2 

4.3 ± 1.8* 

1207 

Gl 

<8 4 

. . . 

1266 

F2 

<4.1 


1280... 

K7 

<2.1 

<4.3 

1284 

A9 

<1.8 

<3.4 

1298 

K2 


<8.0 

1355 

K5 

<3.1 

... 

1380 

A1 


<6.8 

1384 

A4 

13.1 ± 2.1 


1431 

AO 

<2.0 


1454 

K5 

<2.2 


1514 

Gl 

<5.0 


1516 

K6 

<2.6 


1531 

K5 

<5.8 


X-l" 


4 1 ± 10 


(938 939) 

X-2" 

k 

2.8 ± 0.98 


(1053) 

X-3'’ 


5.0 ± 1.5 



* Couni rate Jetermined at the star location. 2 a above local back- 
ground tsee text). 

*’ The X-ray sources not associated with cluster members are at 


a, = 3’‘43'"15’ .5, - 24 01 49" 

sj = 3 43 33 (^, = 24 04 08 

= 3 44 34 6j = 24 13 08 

X-l IS within r of the FO cluster member Hz II 956, but the source 
centroid is also near the field binary system. Hz II 938-939. This system 
is nearer than the Pleiades and contains a K star, so we tentatively 
identify X-l with Hz II 938-939. Hz II 1053 is not a cluster member. 
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Fig 1. — Schematic map of the two IPC fields, showing the location of the surveyed Pleiades stars (points) and the X-ray sources detected (circles). The crosses 
indicate three X-ray sources noi identified with cluster members. Noncluster stars are not shown. 
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Fin 2.—H-R diagram for the Pleiades stars in the two IPC fields. X's represent X-ray deleclions al a level of significance above 2 a. 


Fig. 3.— X-ray detection statistics for the Pleiades stars; «ve show both the 
fraction of the total number of stars in each spectral-type range which were 
detected and the number of individual sources observed in each spectral type 
(given as numbers at the bottom of each bin|. Note that the detection fraction 
shown for the K stars is an upper limit because the Hz II optical catalogue is 
complete only to m, - 14. 


Fig 4.— Summary of X-ray detections, represented by diamonds (2 a level) 
and octagons (3 a levell. For stars observed in both IPC fields, the detections 
(or detection and constraining upper limit) are joined by a solid line. 

Fig. S — Oe)Kndence of the median X-ray luminosity on age for different 
samples of G stars. Bars are located at the values of the median X-ray lumin- 
osity. their le-igth represents the uncertainty in age determination in each 
group The range of observed luminosities is indicated by • (the lower • is 
always fixed by the detection threshold in each group), la) Pre-main-sequence 
G stars (Ku and Chanan 1979, Feigelson and De Campli, 1981), (h) Main- 
sequence G stars in the Pleiades (present survey), (c) Main-sequence G stars in 
the Hyades (Stem et al. 1981 1; (d) Local disk popula ion dwarf G stars (Vaiana 
etal. 1 98 1 , T opka et al. 1 982 ; Rosner el al. 1981). 


FiC. 6 -Maximum-likelihood integral X-ray luminosity functions of 
Pleiades Isolid line} and Hy ades (dashed l-ne) O dwarfs. 


Fig. 7. — Log of the X-ray luminosity vs 5he lug of the projected equatonal 
rotational veloaiy for ihe Pleiades stars later than F7 for which values of. or 
bounds on. these velocity are known Our data are shown as octagons when 
we have measurements for both the values of L, and i sin i. Ik'hen only one of 
the two quantities has been measured Iso that we know an upper or lower limit 
for the other quantity), the corresponding limit is indicated by an arrow point- 
ing in the appropriate direction, when we have only upper lower limit infor- 
mation on both vanables. this is indicated by double arrows. For comparison, 
we also plot (dashed line} the besi fit of Pallavicini et al. (1981) including 
normal dwarf stars. RS O’n's, T Taun stars, and Hyades stars. The solid line 
represents the linear regression (including upper limits, Schmitt 1984) to all 
data without including the RS CVn's logL, s 27.47 + I.48logrsini. 
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ABSTRACT 

We report on the numerical aimulation of two well obaerved SMM 
flares by means of a hydrodynamical code for plasma confined in a 
magnetic flux tube. The parameters of the model loops are obtained 
from the SMM and ground based observations and the only free 
parameters In the models are relative to the spatial and temporal dis- 
tribution of the impulsive energy source. We generally find good 
agreement between the simulated and the observed XRP light curves, 
and we obtain extimatos for the total energy input to the two flares. 
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1. Introduction 

As part of the team E effort in the SMM wrorkshop, we have used the Palermo - 
Harvard Numerical Hydro Model (Peres et al. 1982, Pailavicini et al. 1983, Peres ^ 
OL 1983; see also the chapter on numerical modelling in this book) as a diagnostic 
tool to gam insights into the energetics of the gradual phase in compact loop flares, 
in particular, we have used the model to fit the light curves observed with the X* 
Ray Polychromator on SMM. Since this fitting imposes constraints on the charac- 
teristics of the energy input, we can use the model to evaluate various aspects of 
the flare energetics such as, for instance, the overall energy budget. 

In our approach to the flare modelling, we use the FeXXV light curve, provided 
by the BCS and the PCS instruments, as a diagrrastic toot to ascertain the approxi- 
mate duration and intensity of the impulsive heating term. Our previous work (Pal- 
lavicini et al. 1983) has, in fact, shown that the computed FeXXV line amission 
starts to decline as soon as the heating term is effectively switched off, almost 
independent of the locaMzation of the heat source in the loop. This, however, does 
rwt exclude the possibility of a rapidly-decaying heating component and, in fact, 
the best agreement of the model predictions and the observations has been 
obtained so far by taking into account an exponentially-decaying heating term in 
the decay phase of the FeXXV line. 

To gain information about the localization and the spatial distribution of the 
heating from the SMM ia ages is a more difficult task, in particular, while a bright- 
ening In a pixel (for instance overlaying one of the loop footpoints) is clear evi- 
dence for enhanced emission, it is not necessarily evidence for the localization of 
the heating in that region, unless the dais are resolved in times shorter than the 
conduction travel time. Heat might, in fact, have been deposited elsewhere in the 
magnetic loop and brought there by thermal conduction. For example, tne conduc- 
tion times for the May 7 flare and the Noveaiber 12 flares, which will be discussed 
below, are less than 5 seconds and 20 seconds respectively which should be com- 
pared with FCS observations with a time resolution of, respectively, 150 seconds 
and 70 seconds. Moreover, the coronal part of the loop at the beginning of a flare 
has a relatively low emission measure, and therefore will not emit significantly in 
high temperature lines before the evaporation has filled the loop with large amounts 
of material. 

Both the short heat conduction time and the low initial emission measure in 
corona conspire therefore to give low significance to inferencos on heating 
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localizations cased on simplified linage analysia. More reliable Information can be 
obtained only from images taken erith spatial and temporal resolution high erKMjgh to 
ascertain the loop geometry and to resolve the first few seconds. 

A detailed knowledge of geometry and preflare conditions is crucial for a realis- 
tic flare modelling (e.g., PallavicinI et al. 1083). Comparison with observations has 
been made for the flares of May 7, 1080 at ^ 16:00 U.T. and November 12, 1080 
at ~ 17:00 U.T. These flares have bean chosen because detailed data are avail- 
able both on tfieir geometry and preflare conditions, and because they appear to be 
compact loop flares (our model does rtot consider field line disruption and therefore 
can be applied only to compact loop flares). Ft«r both flares we shall present a 
short summary of the observational data and the results of the comparison with 
model catculntiona. The primary observations were provided by the XRP experiment 
on SMM (Acton et al. 1980). Ha and magnetic data were also u.sed to examine the 
magnetic field structure. 

We have assumed, for simplicity and as a first step In our analysis, symmetric 
transient heating localized at the top of the loop, with a heat deposition having a 
gausslan distribution along the field line coordinate. TMs might be a coarse approx- 
imation In the case of the May 7 flare because Ha observations (Acton et al. 1082) 
indicate that the two foot points had different Intensities at the chromo.spheric level, 
which mig'it be construed as evidence of asymmetric heating in the loop. The 
amount of Impulsive energy deposited In the loop was chosen so as to reproduce 
the uoservcd maximum temperature and the FeXXV Hne maximum flux, while the 
time evolution of the heating was chosen in such a way as to reproduce ^he FeXXV 
line light cuivu. In uidur to achieve high flexibility in changing the model parante- 
ters, we have prescribed the impulsive Ireatktg as a separable function of space 
coordinate and time as 

H(s,t)= H^Xf(t)Xg(s) [1] 


g(8)*exp[ 


2.1 Observations 

The flare on May 7, 1980 at 15:00 U.T. was a compact, short-lived event of 
low total intensity, but a quite energetic one, as indicated by the substantial emis- 
sion in high temperature lines (CaXIX and FeXXV). PCS observations show that it 
was confined within a single 16” x 16" pixel. The light curves were observed by 
the FCS throughout itu entire evolution with low temporal resolution (2.6 minutes; 
see figure 1 for the complete set of FCS observations). An extensive analysis of 
this event, based on simultaneous SMM and Ha observations has been given by 
Acton sLaL (1982), and the X-ray observations have also been discussed by Sim- 
nett (1983). 

As shown by the comparison of Ha and soft X-ray observations, the flare corv 
sis ted of a single compact loop of semilength 6.7 x 10^ cm. The cross section 
of the loop has been derived from the size of the Ha footpoints end is A= 8.8 x 
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10^® cm^, MTtiich implies an aspect ratio (loop diametar/loop length) «■ 0.34. 
Preflare cormal cotiditlons have bean determined using the PCS fluxes observed In 
the aaaie fibcel before the orwet of the flare. A preflare temperature of T« (2.6 ^ 
0.6) X 10^ K has been derived from the ratio of the resonance lines of NeIX and 
MgXI. 


The initial pressure and density have been Inferred from the loop temperature 
and length using the scaling law of Rosner, Tuciter and Valana (1078). The derived 
initM pressure was p^ ~ 10 dyne cm~^, implying a coronal density of n« 1.4 x 
10'^ cm*^ before the flare. This density is near the up, >er bound of the range of 
typical densities observed In active region loops (10* - 10'^ cm'") yet it is cort- 
sistent with the accepted scenario of coronal loops wnich assigns high density to 
small and compact loops and low density to long and extended ones. Cixnparison of 
those Initial values with pressure and density derived from the total em.ssion meas- 
ure in the same PCS pixe* before the flare shows that the latter ones are su'jstarr- 
tiaHy higher than those ^j^rediq^ed on the basis of the scaling law (p > 40 


d\ nos/cm , n « 7x10 


II 


•3 ^ 

cm ); this suggests that the flaring loop w'.« not the 


main contributor to the total emission from this 16“ x 16“ PCS pixel during the 
preflaro stage. This result is not unusual and has been found in many flsres 
ubsorvetl with the PCS; in particulai, it has also been observed in the Nov 12 flare 
discussed b^low. 


2.2 Simulations 

Tlie parameters and the time ovolutioii function used to fit the observe tiorts 
were 

•3 • 1 

H = 1 76 ergs cm sec 
o ft 

s = 5.7 X 10° cm. 

o= 10® cm. 


t/160 

* exp [(150~t)/T] 


O < t < 160 seconds 
t > 160 seconds 


In a first attempt, the decay time r was assumed to be ^ero (i.e., an instan- 
taneous switch-off of the heating). Although our calculation has been performed 
with the transient heating localized at the top of the model loop, we wish to 
emphasize that this choice is not the only possible one. PreMous numerical code 
results (Pallavicini el_flL 1B83) have shown little difference between top and bot- 
tom heated flares, at least as far as predictions of light curves in the various spec- 
tral lines are concerned (see however Pores et al. 1983 for other possible diog- 
nostics). In any case since the May 7 flare was imaged in only a single resolution 
etement of the PCS, the data cannot give any information about the host source 
localization. 

As shown in figure 2, the above set of parameters allows us to tit quite well 
the PeXXV and CaXIX BCS light curves during the riso phase of the flare. Owing to 
the low temporal resolution of the PCS observations of this flaru (2.6 minutes), the 
PCS light curves cannot be used to further constrain the fitting of the rise phase. 
It is apparent that our first cholco of the decay time (r= 0) does not fit the PeXXV 
and CaXIX observations in ttiat the predicted light curves decay faster than 
observed. Pigure 3 shows an analogous diaagreemont between predicted and 
observed OVIli and No IX light curves in the decay phase. Our next attempt has 
been to assume a non- zero decay time for the heating function, all the other param- 
eters remaining the same. Pigure A shows our results for 30 sec. 
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The large Increase in density and in temperature in the course of a flare gen- 
erally causes a decrease of an tiio characteristic time scaies for a given loop 
model (for example, the radiative cooHng time and the conduction cooling time of 
the corona, the sound travel time across the space grid, etc.) and therefore 
imposes small computational time steps. In the case of the May 7 flare this 
requirement eras particularly severe, oxpecialiy during the early decay phase. This 
difficulty and the practical limitations of the available computer time have induced 
us to omdel only the first 20 seconds of the decay phase. It is. iiowever, already 
evident that the heating decay time is still too short, and that a decay time of 
60 seconds would bring our predictions much closer to the observations. 

Among the energy terms considered in our model i.e., thermal, gravitational and 
kinetic, the thermal energy is, os expected, the dominant component by, tipically, 
an order of magnitude throughout tho flarri evolution. Figure 6 shows the evolution 
of the ^otal thermal energy content of the flare during the rise phase. 

According to our model, the total energy deposited in the flare during the rise 
phase is 2.0 x 10^^ ergs. Allowing for an additional energy deposition during the 
decay phase with an exponential docay time of 60 secs., one obtains a total 
energy input of 5 x 10^^ ergs. 

Some uncertainty in the total energy input ia related to our lack of knowledge 
of the site of energy deposition. The energy input necessary for a fit of the PCS 
linos can, in fact, be higher by a factor of, at most, two in the case of energy 
deposition localizad near the loop footpuints. 


3 The No ve mber 12 fla re 


3.1 Observations 

Tho hare of November 12, 1080 at 17:C0 U.T. was substantially ddferent from 
the May 7 event in that It had a larger spatial extent (covering at least three 
1b"x15" PCS pixebi) and a somewhat lower density and pressure in the preflare 
state. Ihe light curves in soft X-ray lines appear to be typical of most flares 
observed with the PCS. Since the minimuPi integration time stop required during the 
decay phase was larger than in the cuso ot the May 7 event, tho simulations could 
be carried further than in tho former case. 

An extensive analysis of this event, based on all thu available ground based 
and SMM observations is being carried on us part of thu UK-SMM work.shop (Mac- 
hidice filAL 1084). In this soction we will refer only to XHP and Har observations as 
f'Mxe pertinent to the gradual phase. 

The magnetic field configuration has been inferred by comparing PCS rasters - 
obtained at 70 secs interval ■ with simultaneous Ho images provided by tiio SOON 
network. The Ha images show loop footpoints on opposite sides of the magnetij 
neutral line, suggesting a loop or a family of loopa arching between the two miiri 
Ha emitting regions. Soft X-Ray images taken by PCS show an elorigatod arc-shaped 
structure, brighter at both ends In the low temperature lines NeiX and MgXi during 
the early phase ot the flare. At the same time, high temperature linns (mainly SXV 
and PeXXV) appear to be stronger at what appears to be the top of a loop. The 
low temporature soft X-Ray regions coincide, within the PCS spatial rnsclution (15**) 
zr.J ine alignment accuracy (10"), with the Ha footpoints. while the PeXXV high 
temperature region is located in between the Ha/NeIX footpoints, showing evidence 
ot enhanced brightening at the loop apex. TNa fact supports o working hypothesis 
of heating localized at the loop apnx. Diagnostic tecniqiies applied to BCS linos 
sht>w that the coronal portion of the loop was initially at a temperature of 2x10' 
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K and decayed rapidly during the flare evolution. 

Froei the combined analyaia of the Ha end PCS observations we And that the 
main flaring structure was sesdiangth L - 2. x 10^ cm and cross- 

sectional area A* 2.6 x 10^' cm^. The aspect ratio is therefore a« 0.14. The 
complex magnetic structure of the region and the spatial distribution of the 
observed Ha emission also suggests that other smaller loops were also probably 
involved in the flare, but that they were not dominant throughout its evolution. 

Preflare conditions were determined using the same method as for the May 7 
flare. The preflare temperature was T= 3.0 ^ 0.6 x 10° K and the preflare pressure 
derived on the basis of the scaling laws of Rosner, Tucker and Valana (1078) was 
p^s 6 dynes cm’^, corresponding to an average density n^- 7x10^ cm'^ in the 
coronal portion of the loop. The density and the pressure derived from the total 
emission mecsure in the same PCS pixels during the preflare phase were substan- 
tially larger than (and not consistent with) the same scaling laws, again suggesting 
that the flaring loop was not the dominant one during the preflare phase and that 
the area imaged into PCS pixels was filled with structures other than the loop 
which later flared. 

The impulsive heating localization and distribution can be inferred from PCS 
observations which show PeXXV emission to be confined within one central 16“x16'' 
pixel, which suggests (aitough this cannot be taken as proof) that the transient 
heating was localized at the loop apex; the size of the pixel where the PeXXV 
emission was recorded at the beginning suggests an upper limit for the dimension 
of the directly heated region and therefore we consider a gaussian spatial distribu- 
tion along the field lines, with 9= 5 x 10° cm. The assumption of symmetrical heat- 
ing appears to be more justified in this case. 

3.2 Simulations 

The light curves for the November 12 flare have been fitted using the same 
procedure as for the May 7 flare. The heating source has the functional form given 
by equation [ 1 ] Mrith: 

H = 10 ergs sec"’ cm'° 

s = 2 X 1^ cm. 
o R 

a= 6 X 10° cm 


f(t) = 


^ exp [ 


( 180- t)/T ] 


0 < t < 180 seconds 


t > 1 80 seconds 


Several cases have been computed for different values of the heating decay 
time (r= 0, 30, 60, 140 secs). The model with r= 60 secs is the one which fits 
the observations best. Notice that the fitting of PeXXV and CaXIX lines is a very 
sensitive function of r, owing to the strong dependence on temperature of the 
emissivity of these lines. 

Fig. 6 shows a comparison of observations and numerical simulation for lines 
observed with the PCS. The observed light curves are obtained summing the contri- 
butions from the three brightest pixels, where a loop-like structure was imaged. The 
observed fluxes in the OVIll line have been multiplied by a factor of two to take 
into account an apparent inflight decrease in sensitivity In this channel (Acton 
1882, private communication). 

As shown in figure 6, fluxes for all the lines, other than NeIX, predicted by our 
model are in quite satisfactory agreement with observations. Not only is the general 
time evolution well-reproduced for all lines other than NeIX, but the absolute flux 
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values also agree with the ohservations within the eatimated accuracy of the FCS 
calibration. The latter is not known, at present, to better than a factor of two (cf., 
Pallavicini et a1. 1083). Notice that the Integrated line fluxes have been obtained 
from our one-dinensional model by aiultlplylno the observed loop cross section by 
the computed loop emission per unit area. No free parameter has been arbitrarely 
introduced to match aiodel predictions with the observatfcma other than the shape, 
the central position and the normalization of the impulsi>'a heating function. How- 
ever the high thermal conductivity in the corona makes the detailed shape of the 
heating function the least important factor among the three. 

The discrepancy between the observed and the predicted fluxes for NeIX is at 
present unexplained. We notice, however, that the ratio of OVIll and NeIX fluxes is 
not consistent with what wniuld be expected from the emissivity functions calcu- 
lated by Mewe et al. and from the assumed chemical abundances (relative 

to hydrogen) of 2.8 x 10'^ and 4.5 x 10* , respectively: tira flux ratios should be 
larger than 10 for all the temperatures greater than T~ 10^ K. Several explana- 
tions are possible, among which a contaminations of the NeIX channel by the 
nearby iron lines is perhaps most likeky. The occurrence of this contamination could 
be investigated in the course of the second SMM mission .by means of spectral 
scans around the line center, in the absence of any better resolution of this 
discrepancy, one could remove the difficulty by arbitrarily decreasing the observed 
NeIX flux by a factor of two; the data presented above, however, have not been 
corrected in this ad hoc way. 

Rgure 7 shows the close agreement between the model predictions and the 
BCS observations of PeXXV and CaXIX resonance line intensities. Notice that the 
BCS observations have been empirically scaled down by a factor 1.6 to bring them 
into agreement with the PCS observations. The BCS, hn fact, integrates over a 
much wider field of view than the PCS pixels where most cf the PeXXV and CaXIX 
is detected. The total amount of energy derived from our siawlation with T= 60 
secs is W= 7.2 x 10^^ ergs. This value, as well as the somewhat smaller value 
derived for the May 7 flaio, falls well within ttie range of observed energies in 
compact flares (Svestka 1981, Simnett 1983). The thermal energy is the dominant 
component of energy, by more tiian one order of magnitude, at any time in the 
course of the flare. In figure 8 we present evolution of the total thermal energy in 
our model. 


4. Discussion and Conclusions 

The above results show that our numerical simulations can satisfactorly repro- 
duce light curves of coronal lines observed in real flares. Por the November 12 flare 
the agreement can be considered good, especially if one takes into account the 
uncertainties due to the chemical abundances, PCS absolute calibrations and 
indeterminacy in the loop cross section. For the May 7 flare, we have been able, so 
far, to fit only the rise phase both of BCS light curves predicting the termal energy 
content and to study the start of the decay phase. Even so, however, the compari- 
son has allowed us to put constraints on the energy deposition process and to find 
agreement with some observations. In particular, we have found that in order to fit 
the flare with our model, additional energy must be provided during the decay 
phase even in this compact event. 

keeping in mind that we use only the PeXXV line intensity light curve as a 
guide to derive the heating evolution and that the modelling depends only on two 
parameters, the transient heating power and its duration. It is interesting to note 
that a good fitting of this line leads directly to a good fitting of all the other FCS 
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and BCS lines. Since the link eaong the Ight curves in our cosyutetion is provided 
by the feet that they are computed from a hydrodynasric model our result ahowr 
that the gradual phase of compact flares Is an essentialy hydrodyiwunlcal 
phenomenon and that the FeXXV earission is actually a tracer of the impulsive heat* 
ing characteristics. While it is obvious that this one-dhaenalonal, one*fluid hydro 
code cannot futty reproduce the complexity of the hare phenoeterwn in all its 
details, these simulations do show that it does provide an overall description which 
has the right macroscopic physics in it and is consistent with the observations of 
the XRP light curves. 

Different considerations apply when we try to consider the influence of small 
scale effects such as, for instance, turbulent inhibition of thermal conduction, rxxt- 
Nnear (but classical) modifications to the Spitzer conductivity, inhomogenelties 
across the field lines and the like. A description of the physics underlying these 
pheixmena, which would allow us to include these effects in our model. Is presently 
lacking; in fact, many of theso basic physical effects are stHI poorly understood for 
laboratory plasmas. Prospects of a full understanding of these aspects for a>stro- 
physical phenomena thus appear at the moment even more remote. 

With these caveats in mind, we can conclude that the numerical simulations 
provided by our code, as well as by other analogous codes, represent a first 
encouraging step in comprehensive flare modelling, and are capable of tackling more 
detailed problems, such as the fitting of line profiles. At the same time, the obser- 
vations from the second SIMM mission should be aimed at providing better diagnos- 
tics for placing more stringent constraints on the basic physical assumptions of 
these models. In particular, high time resolution observations of line intensities and 
line shapes during the first few seconds of a flare should help in identifying the 
site of energy deposition, and should allow a more detailed understanding of the 
dynamic evolution of flare plasmas. 
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Figure cautions 


Figure 1 - PCS observations of the May 7 1980 flare, times o correaponds to 
14:52:00 UX 

Figure 2 - May 7 1 980 flare: synthesi 2 ed light curves (solid lines) from the model 
with an instantaneous switch off of the transient heating at time= 160 seconds, 
and BCS observations (crosses) of Fe XXV and Ca XIX lines. Time = O corresponds 
to 14:54:36 UX 

Figure 3 - Comparison of synthesi 2 ed (solid lines) and observed (crosses) photon 
fluxes for the OVIll and NeIX lines for the decay of the May 7 flare. The model as 
well the time definition are as in figure 2. 

Figure 4 - May 7 1980 flare: synthesized light curves (so'ld lines) from our model, 
in the case of an exponentiul decay of the heating after time= 160 secs with e- 
folding time= 30 seconds, and BCS observations (crosses) for Fe XXV and Ca XIX 
lines. The time definition is the same as in figure 2. 

Figure 5 • Predicted evolution of the total thermal energy for the same model of the 
May 7 flare reported in fiyure 4. 

Figure 8 • Predicted (solid lines) and observed (crosses) fluxes for the six lines 
observed with the PCS for the Nov 12 1980 flare. The model presented here 
corresponds to a continuous heating up to time= 180 seconds followed by an heat- 
ing decaying with an e-folding time of 60 seconds. Time = 0 corresponds to 
1 7:00:00 U.T. 

Figure 7 - Predicted (solid lines) and observed (crosses) fluxes for the Ca XIX and 
Fe XXV lines observed with the BCS for the Nov 12 1980 flare. The model as well 
as the definition of time are the same as in figure 6. 

Figure 8 - Predicted evolution of the thermal energy for our model of the Nov 12 
1980 flare. The model is the same as in figure 6 and 7. 
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